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ALGORITMO PARA A ANÁLISE TOPOLÓGICA DE REDES () 


RESUMO 


São expostos os fundamentos teóricos e a orga- 
nização lógica dum algoritmo destinado à análise auto- 
mática de redes pelo método topológico, usando um 
computador digital. O algoritmo node ser usado indi- 
ferentemente com informação de incidência relativa a 
circuitos ou a cortes, sendo possivel a análise de redes 
com ligação magnética e incluindo elementos não reci- 
procos activos ou passivos. 


1 — INTRODUÇÃO 


À análise de redes lineares, invariantes, de paráã- 
metros concentrados conduz ao estabelecimento de siste- 
mas de equações lienares simultâneas, cuja matriz A dos 
coeficientes é da forma, 


A =:35 PD" (SM) (1) 


em que P é diagonal de dimensão nxn, 5 e D são rec- 
tangulares de dimensão mxn fm< an). 

Os elementos de P descrevem as propriedades fisi- 
cas dos componentes da rede; os elementos de 5 e D 
assumem exclusivamente os valores O, + 1 ou —1 e 
descrevem a forma como esses componentes se encon- 
tram interligados para formar a rede, 

A cada coluna de 5 ou de D e portanto a cada ele- 
mento da diagonal de P, pode fazer-se corresponder um 
elemento dum conjunto de n ligações sobre o qual se 
definem dois grafos à custa das matrizes 5 e D. Para 
tal, as linhas de S e D são interpretadas como vectores 
de incidência orientada dessas n ligações em sistemas 
de circuitos fundamentais, em sistemas de cortes funda- 
mentais ou de cortes nodais, consoante a análise da rede 
tiver sido conduzida pelo método das malhas, dos 
cortes ou dos nós respectivamente [1]. Como é conhe- 
cido [1], as matrizes de incidência deste tipo têm a pro- 
priedade — E definida por Cederbaum [2], (ver abaixo 
em 2. Def. 1) 
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SYNOPSIS 


The theory and logical structure of an algorithm for 
automatic network analysis by digital computer using 
topological methods are described. The algorithm may 
be used indifferentiy with incidence data given for 
circuits or cut-sets. The analysis of networks exhibiting 
magnetic coupling and including active or passive 
non-reciprocal elements is possible. 


O valor das incógnitas dum sistema de equações 
lineares pode ser explicitado em termos dos coeficientes 
do sistema pelo método de Cramer [3] à custa das 
expressões do determinante e complementos algébricos 
dos elementos da matriz À dos coeficientes. 

O teorema de Binet-Cauchy [4] permite exprimir 
o valor do determinante duma matriz da forma (1) em 
termos dos elementos das matrizes factores sem ter de 
efectuar prêviamente o seu produto, utilizando a expres- 
são para det À, 


2 O (ks Kg eso) Pio Does Pim 
(Kyo Ko, co. km) (2) 
sendo, 
titia 165 velo) Gel d(Ky Ko ses di 
- det Dik,, Ko, +.: Ki) (3) 


onde (k,, ks, ... k, ) designa uma das (7) combina- 


ções de m indices de coluna de S e D extraídos dos 
n índices possíveis sem repetição. S ( k,, ks, ...km | e 
D (k, k;, -..k, ) denotam assim matrizes quadradas 
de dimensão mam extraídas de 5 e D referentes às 
colunas designadas pelos índices ky,, ks, -.. k,. O 
produto dos determinantes relativos a um mesmo grupo 
de colnunas de S e D dão assim origem a um coeficiente 
C(ky ko, -. k, |) que em virtude das matrizes S e D 


(*) Comunicação do Laboratório de Medidas Eléctricas do Instituto Superior Técnico sobre trabalhos integrados no Projecto de 
Investigação TLE-5 subsidiado pelo Instituto de Alta Cultura, apresentada no |l Simpósio sobre as Teorias da Informação e dos 
Sistemas que teve lugar na Faculdade de Engenharia da Universidade do Porto em Setembro de 1972. 


(**) D' denota a transposta da matriz D. 
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terem a propriedade — E, só toma os valores O, + 1 ou 
— 1.0 desenvolvimento (2) é portanto uma soma homo- 
génia de produtos p,, P,o.-. Pim Ce m elementos de 
P relativos a cada combinação de indices, afectados 
cada um do correspondente coeficiente dado em (3). 

Note-se que o cálculo do complemento algébrico 


dum elemento 6 pa ca matriz A, 


(— PT der A (4) 


“pa pa 
onde A,, representa a matriz resultante de suprimir 
em À a linha p e a coluna q, reduz-se imediatamente 
ao cálculo do determinante duma matriz da forma (1), 
pois Ná obtém-se simplesmente por eliminação da 
linha p de S e da linha q de D. 

Para gerar a fórmula (2). basta pois: 

a) Identificar as combinações de índices de 
coluna para as quais o coeficiente dado por 
(3) é diferente de zero, 

b) Determinar o sinal algébrico desses coeficien- 
tes. 

A concretização de a) e b) incide únicamente 
sobre informação contida nas matrizes S e D, envolve 
assim únicamente a topologia da rede expressa pelos 
grafos associados a S e D. Em conformidade, designam- 
-se os métodos de análise baseados em fórmulas de 
desenvolvimento dos determinantes do tipo de (2) por 
métodos topológicos, Como é conhecido [1], as com- 
binações referidas acima em a) correspondem 1) a 
árvores, 2) a co-árvores comuns aos dois grafos asso- 
ciados às matrizes S e D, consoante estas matrizes 
representarem as incidências das ligações respectiva- 
mente em 1) cortes, 2) circuitos. 

É extensa a bibliografia àcerca de algoritmos para 
a análise topológica de redes. Alguns desses trabalhos 
estão referidos na bibliografia [5] a [11], outros pode- 
rão ser encontrados nas respectivas referências bibliográ- 
ficas. A maioria dos autores dedica-se exclusivamente à 
resolução do problema ajda identificação e apenas para 
um único grafo! (*) Dos trabalhos referidos, apenas 
Mavyeda et al. [6] e [9] põem o problema com genera- 
lidade para dois grafos distintos, envolvendo portanto 
o problema b) da determinação do sinal, fazem-no po- 
rém sômente para o caso específico da informação ser 
de incidência em cortes. . 

A eficácia dum algoritmo destinado a gerar a fór- 
mula (2) depende da satisfação de certas exigências: 

1) Tem de ser exaustivo, isto é, gerar todas as 
combinações citadas em a). 

2) Não deve gerar duplicados. 

3) Não deve exigir, quer para a eliminação de du- 
plicados, quer por qualquer outra razão, a comparação 
de cada nova combinação gerada com todas as anterio- 
res para decidir da sua aceitação, Mesmo tendo em vista 
a mecanização do algoritmo num moderno computador, 


torna-se proibitivo em capacidade de memória e tempo 
de execução qualquer comparação deste género, Aten- 
da-se a que o número de parcelas da fórmula (2) atinge 
facilmente os milhões.("*) 

Os algoritmos de Mayeda et al. [6] e [9], bem 
como aquele que aqui se descreve satisfazem as referi- 
cas exigências, sendo este último mais simples e aplicá- 
vel sem qualquer modificação a matrizes de incidência de 
qualquer tipo desde que tenham a propriedade — E. 

Em ligação com este problema é importante referir 
os resultados devidos a Bott e Duffin [12] e Duffin 
[13], onde o cálculo de determinantes de matrizes da 
forma (1) é tratado do ponto de vista da álgebra de Gras- 
sman(""*) e se expõem as relações entre esta e a ál- 
gebra de Wang [11]. Encarado sob este ponto de vista, 
o algoritmo que aqui se descreve dá a possibilidade de 
gerar o discriminante não-simétrico (os referidos autores 
consideram apenas o caso simétrico S — D) associado 
ao par de multivectores definidos pelo produto externo 
dos vectores de linha de duas matrizes E como Se D 
que figuram em (1). 


2 TEORIA 


2.1 — Definições 

Def. 1 Matriz — E: Qualquer matriz mxn tal que os 
determinantes de qualquer ordem que dela se 
podem extrair assumem únicamente valores 


O + 10 — 4. 

Nota: Os elementos duma matriz — E (determinantes 
de ordem 1) só podem ter os valores O, + 1 
e —1.. 


Def. 2 Incidência: Uma linha e uma coluna duma ma- 
triz — E são incidentes se o elemento no seu 
cruzamento for diferente de zero. Consoante o 
sinal algébrico desse elemento a incidência é 
positiva ou negativa. 


Def. 3 Redução: Operação que consiste em tomar numa 

matriz — E uma linha e uma coluna incidentes, 
designandas linha e coluna de redução respec- 
tivamente, e somar algêbricamente a linha de 
redução a todas as restantes linhas incidentes 
na coluna de redução por forma a anular todas 
as incidências na coluna de redução excepto a 
da linha de redução. 
Designa-se a matriz resultante da redução por 
matriz reduzida e a submatriz desta, proveniente 
de eliminar a linha e a coluna de redução, por 
matriz restante. 


(*) Nos problemas em que S=D, (redes de elementos com uma única porta),a geração das fórmulas topológicas reduz-se 
apenas ao problema a) da identificação para um único grafo. Não surge o problema b) da determinação do sinal visto que os 


coeficientes (3) são quadrados, portanto sempre positivos. 


(**) O número de árvores dum grafo completo de r» vértices é dado pela expressão v! 2, para + = 10 seriam 100 milhões! 


(***) Também designada álgebra exterior. Veja-se a ref, [3]. 
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Def. 4 Estrutura: Considere-se uma matriz — E rectan- 
gular mxn m<. n. Chama-se estrutura a qualquer 
conjunto de m colunas extraídas das n que for- 
mam uma matriz quadrada mxm com determinante 
não nulo. 


2.2 Teoremas. 


Teorema 1: A matriz reduzida duma matriz —- E é 
ainda uma matriz — E. 


Demonstração: 


A demonstração é conduzida relacionando os deter- 
minantes ca matriz reduzida com os determinantes cor- 
respondentes da matriz inicial que por hipótese só to- 
mam os valores O, + 1 ou —1. 

a) Determinantes extraídos de qualquer conjunto 
de linhas incluindo a linha de redução têm o mesmo 
valor que os correspondentes determinantes da matriz 
inicial, pois diferem quando muito por terem sido 
feitas combinações de linhas que não afectam o valor 
dos determinantes. 

b) Determinantes extraídos de qualquer conjunto 
de linhas que não inclui a linha de redução: 

Se as linhas em causa não foram afectadas na 
redução, é óbvio que os determinantes delas extraídos 
coincidem com os determinantes correspondentes da 
matriz inicial. 

Se algumas das linhas foram modificadas na redu- 
ção, ou a coluna de redução faz parte do determinante 
considerado e então ele é nulo por conter uma coluna 
de zeros, ou a coluna de redução não faz parte do de- 
terminante considerado. Neste último caso acrescente-se 
ao determinante os troços correspondentes da coluna e 
da linha de redução. Aparte permutações de linhas e 
colunas que só afectam o sinal do determinante, pode- 
mos escrever o determinante acrescentado na forma 


+ x x 
O 
O: (5) 
Determinante 
considerado 


Desenvolvendo pela primeira coluna conclui-se que o 
determinante considerado só pode diferir no sinal em 
relação ao determinante acrescentado. Este porém inclui 
a linha de redução e está portanto no caso de a). 

Conclusão: Os determinantes da matriz reduzida 
só podem diferir dos da matriz inicial quanto ao sinal 
ou serem nulos. Satisfazem portanto as exigências da 
propriedade — E. 


Nota: A execução da operação de redução duma ma- 
triz — E não gera para valores dos elementos 
da matriz reduzida outros números que não se- 
jam O, + 1 ou — 1. 


Corolário 1.1: As estruturas da matriz inicial e da ma- 
triz reduzida coincidem e os seus deter- 
minantes têm idênticos valores. 


Demonstração: 


Basta atender a que as estruturas são caracterizadas 
pelo valor dum determinante de ordem m incluindo por- 
tanto todas as linhas e que caem por isso no caso a) da 
demonstração do teorema. 


Nota: A execução de operações de redução sobre uma 
matriz não altera a informação nela contida no que 
respeita às estruturas. 


Teorema 2: Designando: 


E, (Mi o conjunto de estruturas da matriz 
; M que incluem a coluna k, . 
E (Mi 5 o conjunto de estruturas da matriz 
* restante da redução de M pela 
coluna k, utilizando a linha I, nela 
incidente. 


O teorema afirma que, 
EM [hj >< EM) (*) (6) 
k: 


Ê 


Demonstração: 


Qualquer determinante de ordem m — 1 da matriz 
restante pode ser acrescentado com a coluna e linha 
de redução e levado à forma (5) quando muito com uma 
mudança de sinal. Ser o determinante considerado dife- 
rente de zero implica que o determinante acrscentado 
também seja e reciprocamente. Em virtude do corolário 
1.1 o determinante acrescentando, que é um determinante 
de ordem m da matriz reduzida coincide com o deter- 
minante correspondente da matriz inicial. Conclui-se 
assim que qualquer estrutura da matriz restante acres- 
centada a coluna k, dá uma estrutura da matriz inicial 
e riciprocamente a qualquer estrutura da matriz inicial 
contendo k, corresponde a existência duma estrutura na 
matriz restante formada pelas restantes colunas. 


Corolário 2.1: Entre os determinantes D duma estrutura 
da matriz inicial e D' da estrutura corres- 
podente da matriz restante existe quando 
muito uma diferença de sinal expressa 
por, 


B=(— pt+Ere py (7) 


(*) O sinal «x» denota a composição dos elementos do conjunto à esquerda (neste caso um único) com os elementos do conjunto 


à direita. 
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em que: 


e Toma o valor O ou 1 consoante a inci- 
dência da linha na coluna de redução 
+ Qu—. 

! Denota o número de ordem da linha de 
redução no determinante D. 

c Denota o número de ordem da coluna 
de redução no determinante D. 


Demonstração: 


Basta considerar o determinante D desenvolvido 
na matriz reduzida pela coluna de redução que tem um 
único elemento não nulo (incidência na linha de redu- 
ção). O menor relativo a este elemento é precisamente 
o determinante D'. 


Nota: Os números e e só dependem da linha e colu- 
na de redução consideradas e são por isso os 
mesmos para todas as estruturas da matriz res- 
tante. Sômente c poderá ser diferente de estru- 
tura para estrutura em virtude de depender da 
posição que irá ocupar, no determinante D corres- 
pondente, a coluna de redução. 


Teorema 3: Seja |, uma linha de Mek,, ks, .. k,0 


conjunto das colunas incidentes nessa linha. 
Designando por M,. 4 k a matriz 
I o, “ . . p 


que resulta de omitir as colunas designadas 
em indice na matriz M, tem-se que, 


Ei M) Ex, UM) + Eu, (My ,) e. Es (MA, ksccksp) 


(8) 


onde os conjuntos do segundo membro 
são disjuntos. 


Demonstração: 


Os conjuntos que figuram no segundo membro 
são ôbviamente sub-conjuntos do conjunto do primeiro 
membro e são disjuntos uma vez que em cada parcela 
se omitem as colunas que figuram obrigatóriamente nas 
parcelas anteriores. Basta então provar que não existe 
nenhum elemento de E(M) que não esteja incluído num 
dos conjuntos do segundo membro. 

Ora um elemento de E(M) que não estivesse 
incluído no segundo membro não conteria nenhuma das 
colunas k,, ks. .. k, incidentes na linha |, e o seu 
determinante seria nulo por conter uma linha de zeros. 
Esta conclusão contradiz a hipótese de se tratar dum 
elemento de E(M) e portanto uma estrutura. 


Corolário 3.1: A identificação das estruturas duma ma- 
triz M de m linhas pode farzer-se à custa 
da identificação das estruturas dum certo 
número de matrizes de m-—1 linhas deri- 
vadas da primeira por redução sem gerar 
duplicados. 


Demonstração: 


Basta combinar os resultados expressos pelos teo- 
remas 2. e 3. À parcela genérica do desenvolvimento (8) 
pode escrever-se utilizando (7), 


/, 
Ee Mp lpeeodo r (6, )> EM bpeecbo Pk Do cas 


Como estas parcelas são disjuntas, a identificação 
por esta via não dá origem a repetições. Note-se que a 
linha de redução , é a mesma para todas as parcelas, 
representando k (p |, ... Ss) O conjunto de colunas 
incidentes nessa linha. 


Os teoremas que se seguem permitem generalizar 
os resultados anteriores para o caso em que estão 
em causa duas matrizes S e D da mesma dimensão 
mxn tm < in). 


Teorema 4: Seja k,, ks. ... k o conjunto de colunas 
incidentes numa determinada linha | j da ma- 
triz S ou da matriz D . Tem-se que, 


ES MN ED) E, (SIM Es, tD) (9) 
Ea, Se) NM Eç,iD,) 
FER Chohochaca) NE Depkdckis! 
Demonstração: 
Designando por k,, ks, ... k, um conjunto de 


q=>n-m+ It colunas das matrizes S e D, podemos 
decompor os conjuntos de estruturas de S e D em 
parcelas da forma, 


E1$) Exp 1S1 + EglSk,) A TER Skoda) 
ED) E, (D) E, »1Dy 1 praca E, (De, ko...k 


Note-se que à a excepção das parcelas dum e doutro 
cesenvolvimento referentes aos mesmos indices, todos 
os conjuntos parcelares dos segundos membros são 
Cisjuntos. Fazendo a intersecção membro a membro 
obtém-se simplesmente, 


EIS) MN E(O) Ex, (SIN E, 1D) 


RS DRE PL E 

Supondo agora que as colunas de índices k,, ks. 

- kk, correspondem exactamente às incidências 
numa certa linha, por exemplo da matriz S, teremos 
que todas as parcelas da decomposição de E(S) além 
das que correspondem a estas colunas serão vazias 
(Teorema 3). Basta portanto na expressão da intersecção 
considerar apenas as parcelas relativas a esta mesma 


colecção de indices. 
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Corolário 4.1: A identificação das estruturas comuns 
a um par de matrizes S e D da mesma 
dimensão com m linhas pode fazer-se à 
custa da identificação das estruturas 
comuns dum certo número de pares de 
matrizes com m — 1 linhas, derivadas de 
Ss e D por redução duma certa coluna, 
sem gerar duplicados. 


Demonstração: 
Procedendo do mesmo modo que se procedeu para 


o corolário 3.1, obtém-se para parcela genérica do 
cesenvolvimento (9), 


E SK, acto! PM ER Da, ok 


| 
pl 


e | mt; M Ei | OR 
ad ES ky kpeka] da Dept puckpo 


Nota: Ao levar a cabo as reduções necessárias para a 
formação de cada uma destas parecelas, a coluna 
de redução k, é a mesma para ambas as matrizes, 
& linha de redução para uma das matrizes 
é a mesma em todas as parcelas l, para a outra 
pode ser qualquer linha incidente na coluna de 
redução. 


A propósito das expressões de decomposição dos 
teoremas 3 e 4 não foi explicitamente considerada a 
possibilidade do número s de colunas incidentes na 
linha L, poder ser superior a n— m + ?. Quando assim 
acontecer, as parcelas relativas às colunas de índice 
k comp -n-—m+ 7 são evidentemente vazias, pois 
que feita a supressão de colunas já não restam suficien- 
tes para formar uma matriz quadrada m x m. 

Finalmente atenda-se ao problema b) da determi- 
nação do sinal dos coeficientes C de (3). Utilizando o 
resultado do corolário 2,1 tem-se, 


Os números c, e c, que se referem à ordem da 
coluna de redução numa e noutra matriz são iguais 
visto que a coluna de redução é, como se salientou 
acima, a mesma nas duas matrizes. O sinal algébrico 
de C não é por isso afectado por estes números, podendo 
resumir-se o resultado anterior a, 


À E E, T+ É, , | | 
C-Ds. Do=[(—1) HD. Do 


(11) 
Conclui-se assim que o sinal algébrico a introdu- 


zir para exprimir o produto de determinantes de ordem 
m em termos do produto dos determinantes de ordem 
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m-— 1 é o mesmo para todas as estruturas duma certa 
parcela (10) do desenvolivmento (9). A inserção deste 
sinal pode assim ser feita automáticamente associan- 
do-o ao elemente k, que é posto em evidência em 
cada parcela (10). 


3 ALGORITMO 


3.1 — Descrição 


A aplicação sucessiva da decomposição expressa 
pelo teorema 3, (teorema 4 para o caso de duas matri- 
zes), permite reduzir sucessivamente o número de linhas 
das matrizes consideradas até se atingir uma única 
linha, As estruturas neste caso são simplesmente forma- 
das pelas colunas incidentes nessa última linha, (colu- 
nas incidentes simultâneamente na linha restante duma 
e doutra matriz no caso de duas matrizes). 

A identificação e determinação dos sinais das es- 
truturas pode assim fazer-se através dum algoritmo 
recursivo capaz de efectuar as operações de redução 
relativas à transição do nível m para o nivel m — T. 
O encadeamento dos sucessivos níveis pode ser ilustrado 
por um diagrama em forma de árvore como o da Fig. 1 
que ilustra a redução duma matriz de cortes nodais do 
grafo representado na mesma figura, A matriz inicial 
com três linhas dá origem por aplicação do teorema 3 
a uma descendência de matrizes de duas linhas e cada 
uma destas finalmente a uma descendência de matrizes 
com uma única linha. Os elementos não nulos destas 
últimas dão os nós terminais do diagrama. Na figura 
cada seta a cheio denota assim uma operação de redu- 
ção da coluna assinalada no nó terminal respectivo, 
Junto do Ínicio de cada seta estão marcadas as colunas 
suprimidas antes da redução e ao longo da seta a 
linha de redução usada na mudança de nivel. 

À cada nó terminal do diagrama da Fig. 1 corres- 
ponde uma estrutura que é formada pela coluna corres- 
pondente ao nó terminal mais todas as colunas de re- 
dução assinaladas nos nós ao longo da sua linha de 
descendência a partir do nó superior. 

As setas a tracejado indicam o progresso do algo- 
ritmo através do diagrama. Esse progresso faz-se sis- 
temâticamente de cima para baixo e da esquerda para 
a direita. O progresso em profundidade tem prioridade 
a só termina quando se atinge o nível mais baixo, 
seguindo-se então a exploração horizontal a esse nível. 
Exausto um nivel o algoritmo transita para o nível ime- 
diatamente acima. Caso não esteja ainda exausto prosse- 
que a exploração horizontal tomando a coluna de redu- 
ção imediatamente a seguir à última que foi tomada a 
esse nível e passa de novo aos níveis mais baixos afim 
de explorar toda a descendência dessa operação de redu- 
ção. O algoritmo termina quando, exaustos todos os 
níveis, ocorre o regresso ao nivel zero, 

As linhas e colunas usadas nas reduções têm de ser 
marcadas de modo a não serem usadas a niveis inferio- 
res áquele em que foram tomadas. Quando se dá a 
transição para um nivel superior essas marcas são apa- 
gadas no que se refere a niveis inferiores ao corrente. 
Deve nortar-se que sempre que se efectua uma combi- 
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k1 ka ka k4 ks 
11 1 0 1 1 
lo |-1 0 0 
13 0 -1 Õ -1 
DIAGRAMA DL REDUÇÃO 
(n) k1,k3 i 
14 dá k1 ; 
/ 1 
5 
(1) 
Z ho rka3 
7 ka 
“A 
do do dA Es 
Ad 
ONBRONRO 


Expressão escrita pelo algoritmo: 


El EkoCEks kg) tka( ks 


Soma homogénia de produtos: 


kk (ER) kaka +ks ERG) + ks + ko (+ kg 


Fig. 1 


nação de linhas para efeitos de redução se tem de fazer 
a soma a toda a extensão das linhas e não apenas sobre 
as colunas ainda livres. Deste modo preserva-se sempre 
a informação contida na matriz relativa a estruturas 
e as sucessivas transformações podem ser sobrepos- 
tas sem necessidade de preservar a forma anterior da 
matriz. Esta circunstância não só simplifica a programa- 
ção do algoritmo como o torna pouco exigente em capa- 
cidade de memória. 

Uma das características mais interessantes deste 
algoritmo está na estreita ligação que existe entre a se- 
quência de operações que realiza e a sequência dos sim- 
bolos duma fórmula algébrica para o valor de det A, 
veja-se (1). Com efeito, pode obter-se essa fórmula 
fazendo corresponder, à transição para um nível mais 
baixo a abertura dum parêntesis, à redução da coluna k, 
a inserção do elemento + p; correspondente da matriz 
P acompanhado do sinal algébrico conveniente calculado 
de acordo com (11) e finalmente à subida de nível 
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— Exemplo demonstrativo do funcionamento do algoritmo 


o fechar dum parêntesis. Na Fig. 1 pode ver-se a ex- 
pressão tal como seria «escrita» pelo algoritmo para o 
exemplo dado, (os sinais são todos positivos por se 
tratar duma única matriz, isto é, S= D na expressão 
(1). Compare-se esta escrita com a escrita da mesma 
expressão na forma duma soma homogénia de produtos, 
cada um relativo a uma estrutura, e verifique-se como o 
algoritmo põe automáticamente termos em evidência re- 
Cuzindo o número total de multiplicações, (efectiva- 
mente 8 contra 16)! 


3.2 — Organização lógica. 
O diagrama da Fig. 2 e a respectiva legenda ilustram 
a organização lógica do programa de execução do algo- 
ritmo. 
Na «selecção de 
de escolher, 


linha» empregou-se o critério 
entre as linhas ainda livres de S e D, 
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Princípio 


(1) | Inicialização do 
indicador de nivel 
N=0 


Wao 


Selecção da linha Inscrição de:| 
de 
reduçao 


(10) 
Subida de Selecção da coluna e| pescida de 
nivel = sua | | nivel 
N = N-1 Nao redução  N+1 


extate 


(11) 
Marcação das 


linhas 


DPesmarcação 
de linhas e 
colunas 


Inscrição de: 
) 


Fim 


Bloco [2): Selecção, entre as linhas livres (não marcadas) de ambas as matrizes de incidência, de uma 
linha que tenha o menor número de colunas incidentes. À linha seleccionada fica associada, 
como linha re redução, à transição do nivel corrente para o imediato. 

Caso o número minimo seja zero, (não existe), significa que o ramo de descendência corrente 
estã vazio, devendo transitar-se ao nivel anterior afim de explorar outra linha de descendência. 


Bloco (3): Selecção da primeira coluna livre (não marcada) incidente na linha de redução corrente. 
Busca duma linha incidente nesta mesma coluna na outra matriz, Redução de ambas as ma- 
trizes e determinação do sinal, expressão (11). 
Caso já não existam colunas susceptíveis de redução significa que o nivel está exausto, 
devendo transferir-se para o nível anterior. 


Bloco [4] e (7): Sempre que uma coluna é tomada para redução é marcada (bloco 4). Antes de transitar 
a um nivel mais baixo, as linhas de redução do nivel corrente são igualmente marcadas. 
Às marcas são tais que permitam identificar o nivel a que foram estabelecidas. 

Bloco (12): Apagamento das marcas relativas a níveis inferiores. 

Blocos (5), (9) e (13): Os símbolos apropriados são acrescentados à expressão da fórmula, ou são 
efectuadas as operações correspondentes sobre números ou expressões representivas dos 


parâmetros Pj - 


Fio. 2? — Organigrama do algoritmo 
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aquela que tem o menor número de incidências. Isto é, 
aquela que dá origem à descendência menos numerosa 
no nível imediato. Este proceder aumenta a eficiência 
do algoritmo mas não é essencial, Qualquer linha 
livre pode ser usada como linha de redução para o esta- 
belecimento do nível seguinte. 

Tal como está ilustrado na Fig. 2, o algoritmo 
«escrever a fórmula do determinante — blocos 5, 9 e 
13 —, normalmente estes blocos seriam substituídos 
por outros em que se efectuassem as correspondentes 
operações numéricas ou algébricas sobre os elementos 
da matriz P. Neste caso a informação produzida pelo 
algoritmo não necessita ser retida em memória e é 
utilizada à medida que é produzida. 

Excepto em problemas muito simples, a impressão 
da fórmula não tem interesse por ser muito longa, con- 
vém porem ter presente a possibilidade de a armazenar 
sequencialmente em banda magnética para realização 
subsequente dos cálculos utilizando a informação regis- 
tada na mesma sequência. Uma vez estabelecida a fór- 
mula relativa a uma certa topologia da rede ela pode ser 
usada para cálculos relativos a quaisquer valores dos 
parâmetros da rede sem necessidade de repetir a análise 
topológica. 


4 — CONCLUSÕES 


Estabeleceu-se a teoria e o programa de mecaniza- 
ção dum algoritmo simples e eficiente para a obtenção 
do desenvolvimento de determiantes não simétricos que 
ocorrem na análise de redes. 

O algoritmo opera simbólicamente com base na 
informação topológica traduzida por um par de mastri- 
zes — E de incidência estreitamente ligadas ao esquema 
eléctrico da rede analisada. O algoritmo é aplicável sem 
modificação quer a informação topológica diga respeito 
a circuitos quer a cortes 

Evidenciou-se uma correspondência completa entre 
os sucessivos passos do algoritmo e a sequência de 
símbolos que traduzem uma fórmula algébrica capaz de 
fornecer o valor do determinante em função dos valores 
dos parâmetros da rede. É digno de nota o facto do al- 
goritmo automáticamente colocar termos em evidência 
reduzindo o número total de operações de multiplicação 
a efectuar ao realizar os cálculos. 

O algoritmo, que é fundamentalmente um processo 
de gerar uma fórmula, permite a utilização imediata da 
informação gerada para a realização dos cálculos, sem 
necessidade de armazenamento dessa informação e per- 
mite igualmente o seu registo sequencial e utilização 
posterior. 


Com base no organigrama da Fig. 2 foi escrito 
um programa em FORTRAN IV que foi ensaiado com 
sucesso num computador digital.(') 
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NOS . 
LABORATÓRIOS 
SUBTERRÂNEOS DA 
SKF UM CIENTISTA OCU- 
POU-SE DURANTE ANOS A APER- 
FEIÇOAR A MÁQUINA DESTINADA A RE- 
VOLUCIONAR OS FUNDAMENTOS DA TECNOLOGIA 


DOS ROLAMENTOS. 


Se não parecemos preocupados, devíamos parecer. Este aparelho vai causar 
grande agitação entre aqueles que lidam com rolamentos e são muitas 
mais pessoas do que as que podemos imaginar. Este é o primeiro 
sistema que permite o controle do estado dum rolímento e que 
pode prever com certa antecedência uma futura avaria. Em 
primeiro lugar, vai ajudar o utilizador dos rolamentos, 
pois irá indicar quando e porque se devem substi- 
tuir. Vai ajudar, tambem, os nossos concor- 
rentes, pois os seus clientes podem avi- 
sa-los com uma certa antecedência 
de que vão ter que substituir 
um rolamento. Enfim, ama- 
nutenção preventiva a- 
tingiu uma nova era. 

Se quer saber 
mais porme- 
nores, in- 
forme- 

-se na 
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PRAÇA DA ALEGRIA, 66-A TELEF.: 362301 TELEGAR.: ESKAEF RUA DELFIM FERREIRA, 604 TELEF,: 092054 TELEGR:: ESKAEF- 


TECNICA XVHI 


Técnica 417 Novembro de 1972 


Ano XLMVII 


Volume XXXV 


C. D. U. 691.327: 620.197.14 


PROTECÇÃO DE DEPÓSITOS E OUTRAS ESTRUTURAS DE BETÃO 


POR MEIO DE REVESTIMENTOS 


RESUMO 


Devido ao uso cada vez maior do betão de cimento 
portland tanto no dominio da construção e obras publi- 
cas como no das técnicas e operações industriais, in- 
cluindo a armazenagem de produtos quimicos e alimen- 
tares os mais diversos, estão continuamente surgindo 
novos problemas do seu ataque químico, ao mesmo 
tempo que vão aparecendo no mercado novos produtos 
para a sua protecção cuja adaptação a dadas condições 
de corrosão deve ser antecipadamente confirmada em 
cada caso prático específico, 

No presente trabalho bibliográfico fazem-se algumas 
considerações gerais sobre o ataque quimico do cimento 
por produtos de uso industrial e afloram-se os processos 
mais viáveis para proteger as superficies de betão por 
meio de revestimentos anticorrosivos. No capitulo final 
resume-se a matéria delineada e indica-se a principal 
bibliografia a que, em principio, se poderá recorrer para 
definir e estabelecer o processo, ou processos mais 
aconselháveis, para revestir o betão de cimento portland 
em face de cada produto industrial ou alimentar. 


1 — INTRODUÇÃO 


Há já algumas décadas que o betão de cimento 
portland tem encontrado aplicações cada vez mais nume- 
rosas na construção de grandes recipientes — reserva- 
tórios, silos, cuvas, depósitos diversos — destinados 
a armazenar não só a água potável ou não potável como 
também diversos resíduos industriais e muitas substân- 
cias de natureza química ou alimentar. Mas grande parte 
destes produtos tem uma acção corrosiva mais ou me- 
nos intensa sobre os constituintes do betão, em espe- 
cial sobre o cimento hidratado, acção esta cujo resultado 
final pode ser a decomposição do ligante hidráulico com 
a consequente alteração da superfície do betão, ou a sua 
expansão e desagregação a prazo mais ou menos curto, 
Ainda que os componentes do betão — cimento, areia, 
agregados — não possuam toxicidade apreciável, o 
possível ataque pelo produto armazenado provocará a 
alteração deste e a consequente modificação da sua 
composição, aspecto, olfacto ou gosto. Esta alteração 
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ACÁCIO DA CONCEIÇÃO | 
Eng. Quimico, Especialista do ST : 
Serviço de Materiais de Construção 


SYNOPSIS cg 


Owing to the increasing use of portland cement 
concrete not only in buildings and public works but in 
process engineering operations — including the storage 
of chemical and food products —, new problems are 
constantly arising in connexion with chemical attacks 
on these structures. At the same time new products for 
protection against these attacks are constantly appearing 
in the market and their efficiency in given corrosion 
problems must be confirmed beforehand in each 
particular case. 

The present bibliographical report begins with 
some general considerations on the chemical attack 
of cement by industrial products, and briefly introduces 
the most feasible methods of protecting concrete sur- 
faces by means of corrosion-resistant coatings. The 
final chapter sums up the subject of the report, and 
lists main references which in principle should be used 
tor choosing the most suitable coating procedures for 
protection of portland-cement concrete against given 
industrial and food products. 


pode ser muito significativa e assaz nociva para o pró- 
prio material armazenado, especialmente se se tratar 
de produtos químicos comercialmente puros e de substân- 
cias de natureza alimentar. 

Se nos recipientes metálicos há o problema essen- 
cial da protecção contra a corrosão ferrosa, nos depó- 
sitos de betão temos que salvaguardar igualmente o 
ataque químico do ligante hidráulico pelo produto arma- 
zenado, com a agravante de que, se for um líquido, ele 
poder penetrar, mais ou menos fãcilmente, através das 
paredes do depósito devido à porosidade e permeabi- 
lidade do betão. As técnicas de protecção interior dos 
depósitos de ferro e aço por esmaltagem cerâmica, ou 
com revestimentos de metais mais nobres que o ferro — 
cromagem e niquelagem, anodização, metalização por 
projecção, pinturas galvânicas a zinco, etc. — não são, 
evidentemente, praticáveis nas superfícies de betão; 
mas não quer dizer que não haja meios viáveis para as 
proteger contra qualquer agressividade química, com 
relativa segurança e continuidade, por meio de revesti- 
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mentos adequados a cada caso prático. No entanto, 
para uma armazenagem prolongada, a temperaturas 
muito fora do normal, não se pode talvez assegurar uma 
absoluta inoquidade a certos produtos químicos ou ali- 
mentares, tais como extractos em fermentação, certos 
ácidos inorgânicos e bases bastante concentradas, etc. 

No capítulo 4 indicam-se, a traços gerais, as dife- 
rentes técnicas para proteger o betão contra o ataque 
químico, de acordo com a bibliografia mais recente 
sobre o assunto. Antes porém, referem-se sucintamente 
os mecanismos gerais da corrosão do betão por produtos 
de uso industrial como princípios básicos indispensáveis 
a uma melhor apreensão do assunto. 


2 — CONSIDERAÇÕES GERAIS SOBRE O ATAQUE 
QUÍMICO DO BETÃO 


Sabe-e que o betão de cimento é um material rela- 
tivamente pouco nobre no que respeita à estabilidade 
a produtos e reagentes químicos, embora possua resis- 
tência mecânica apreciável, incombustibilidade, e suporte 
razoavelmente a acção da água normal e salgada e de 
certos reagentes de alcalinidade acentuada. As causas da 
destruição do betão são numerosas e variadas. Em pri- 
meiro lugar há uma série delas que se relacionam com 
a preparação e o tratamento sofrido durante os primeiros 
períodos a seguir à aplicação em obra — dissecação 
muito rápida, gelamento, dosagem muito rica ou muito 
pobre em cimento, carência ou excesso de fluidez, com- 
posição granulométrica defeituosa, fraca estabilidade 
físico - química dos agregados, etc. —, causas estas que 
podem provocar uma excessiva heterogeneidade e fissu- 
ração do betão, isto é, torná-lo mais vulnerável aos 
agentes corrosivos e contribuir, em grande parte, para o 
seu ataque e desagregação. 

Por outro lado, a destruição do betão é também 
motivada pela fraca resistência química dos compostos 
integrantes do cimento hidratado; instabilidade esta 
que se atribui mais geralmente à cal, à magnésia ou ao 
gesso mas se amplia facilmente aos aluminatos e silicatos 
hidratados, que se podem decompor pela acção de 
certos agentes agressivos em circunstâncias e ambientes 
adversos. 

Esquemáticamente, a hidratação do cimento portland 
pela água da amassadura do betão pode traduzir-se 
quimicamente como segue: 

+ nH,0 


-» 


1) CS + Ca0O livre C,s. n H,0+ 2 Ca (0H) » 


+ 2 Ca (0H), ———. 
ária C, A.6H,0 + C4 Fim' H,0 


C, A.6H,0 C,A.350,Ca, 3] Hs (0) 
Ca F.m H,40 C;F.3S0,Ca. m H, 0 
Sulfoaluminato e sul- 
foferrito de cálcio, 
hidratados 
Cas 
3) C.s +p H,0—— C,S.p' H,0 + Ca (0H), 
: Silicatos Hidróxido 
hidratados de cálcio 
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A cal livre e algum silicato tricálcico são primeira- 
mente hidratados (R. 1), indo o respectivo hidróxido de 
cálcio, libertado na reacção, regular a hidratação do 
aluminato tricálcio (C, A) e do aluminoferrito tetracál- 
cico (C,AF) e a sua combinação com o gesso de presa 
para formar os respectivos sulfoaluminato e sulfoferrito 
de cálcio, hidratados (R. 2). Na continuação da hidra- 
tação e endurecimento do cimento, o silicato tricálcio 
(C,S) e o silicato bicálcico (C,S) são transformados 
progressivamente em silicatos hidratados, com libertação 
do hidróxido de cálcio que torna o meio do betão fran- 
camente alcalino (pH = 12,4) onde esses silicatos 
cristalizam e são estáveis. Mas este hidróxido de cálcio 
(leite de cal) — que continua a libertar-se até final do 
endurecimento — é um composto muito vulnerável ao 
ataque químico e ao arrastamento para fora do betão, 
o que pode motivar o abaixamento do pH do meio, com 
a consequente instabilidade dos silicatos e aluminatos 
hidratados nele inseridos. Nestas circunstâncias serão 
cecompostos facilmente pelos agentes agressivos, es- 
pecialmente os de natureza francamente ácida e certos 
sais e sulfatos alcalinos e alcalino-terrosos em concen- 
tração acentuada. 


De modo que entre as substâncias mais nocivas 
para o betão de cimento podem citar-se: (a) os ácidos, 
mesmo fracos, especialmente os inorgânicos; (b) os 
compostos de magnésio solúveis, especialmente o sul- 
fato; (c) outros sulfatos e cloretos tais como os alca- 
linos e os de cálcio, mas em soluções concentradas; (d) 
as águas puras e agressivas para a cal (CO, agressivo), 
ou pelo, contrário, as muito salinas tal como as sulfata- 
das ou selenitosas; (e) certos compostos amoniacais e 
alguns nitratos; (f) alguns ácidos gordos de óleos ve- 
getais e animais, especialmente os líquidos, tal como 
o oleico e outros ácidos gordos insaturados; (g) alguns 
óleos minerais tal como os «fuel-oils»; (h) certos pro- 
dutos químicos tais como os hipossulfitos e bicromatos 
de potássio; etc. 


No entanto há um certo número de materiais que 
não atacam notóriamente o betão, entre os quais se 
podem incluir os carbonatos e a maior parte dos nitra- 
tos, alguns cloretos e fluoretos, e os silicatos. Também 
a água de cal e outras soluções alcalinas não são agres- 
sivas para o betão desde que não excessivamente concen- 
tradas. Os produtos derivados do petróleo, quando 
isentos de aditivos acídicos, e alguns destilados do 
alcatrão de carvão, não têm efeito muito apreciável sobre 
o betão compacto e bem endurecido, além de certa des- 
coloração, que pode ou não ser indesejável. 

A deterioração do betão manifesta-se, geralmente, 
pelo aparecimento de fendas, primeiro muito finas e qua- 
se imperceptíveis, mas que aumentam progressivamente, 
Por aí penetram mais fâcilmente os agentes agressivos 
que se combinam rápidamente com o hidróxido de cál- 
cio libertado na hidratação ou provocam o seu deslava- 
mento e o abaixamento do pH do meio. O ataque do 
aluminato e dos silicatos de cálcio, hidratados, não 
se fará esperar, especialmente se o agente agressivo 
for apreciavelmente acídico. Se os produtos da reacção 
forem solúveis ou emulsionáveis no meio corrosivo serão 
facilmente exsudados e o desequilíbrio da reacção pro- 
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voca aceleração da corrosão, Como se compreende, 


dos. No entanto, se o cimento for muito de as, EMA 
o estado fisico do agente corrosivo para o betão tem 


sulfatos alcalinos podem também reagir com o aluminato Pa 


4Ca0O. Al0 4. Aq + 2 Na (0H) > 


também um significado muito apreciável devido à capa- 
cidade de penetração e de reacção. Os sólidos secos não 
atacam normalmente o betão também seco, mas alguns 
atacam o betão húmido. Um produto sólido, mas húmido, 
pode atacar o betão, Os gases secos, se forem agressi- 
vos, podem roubar suficiente humidade ao betão para 
tornar posível o seu ataque químico, mas serão bem 
mais destrutivos se estiverem húmidos. Os agentes cor- 
rosivos no estado liquido penetram e reagem muito mais 
intensamente do que se forem sólidos, etc. 

Muitos agentes corrosivos para o betão alteram 
a sua composição química por mecanismos que são 
apenas parcialmente conhecidos ou incompletamente 
compreendidos. A própria água do mar, talvez devido à 
presença de sulfatos, pode ser corrosiva para os compos- 
tos hidratados do cimento, mas destrói mais intensa- 
mente os betões permeáveis ou os que são fabricados 
com cimentos de alto teor em aluminato tricálcico. Por 
outro lado, alguns compostos orgânicos constituídos por 
poli-oxidrilos (glicol, gliceróis, glucoses, ácidos orgânicos 
poli-hidroxilados, etc.) podem atacar lentamente o betão, 
mais provâvelmente através da sua combinação com o 
hidróxido de cálcio libertado na hidratação do cimento. 
A título de exemplo, e muito rápidamente, indicam-se a 
seguir algumas reacções elucidativas de ambos os me- 
canismos principais de agressividade para o betão — a 
acção sobre o hidróxido de cálcio libertado por hidrata- 
ção do cimento (Ca (0H),) e as reacções, expansivas 
ou não, sobre os próprios aluminatos e silicatos hi- 
dratados que o constituem. 

— Acção dos ácidos: Conquanto alguns ácidos 
muito fracos não tenham acção notória sobre o betão 
bem endurecido, todos os ácidos fortes, especialmente 
os inorgânicos, decompõôem os silicatos e aluminatos 
hidratados do cimento endurecido em sílica e alumina 
gelatinosas e sais de cálcio mais ou menos solúveis na 
água: 


3 CaO. 2 SiO,. Aq+ , HCl-——> 3 Ca Cl; + SiO,. aq 


3 CaO. Al,0,. Aq + 5H, 50, — q€a S0,+ Als05 

— Mcção das bases: Em solução relativamente 
diluída (até cerca de 10%) as bases pouco ou nada 
atacam o betão que possui, por si, basicidade acentuada. 
O componente integrante do cimento mais débil a uma 
acção alcalina — a aluminato tricálcico — está normal- 
mente em proporção mais fraca e reage parcialmente 
com o gesso durante a presa do cimento, Mesmo assim, 
o cimento de tipo aluminoso pode reagir com os hidróxi- 
dos alcalinos, por exemplo, formando aluminatos alcali- 
nos bastante solúveis: 


— Mcção dos sais: Manifesta-se geralmente so- 
bre os componentes hidratados do cimento, mas alguns 
deles, especialmente os sais alcalinos, actuam fãácil- 
mente sobre o hidróxido de cálcio libertado na hidrata- 
ção, de preferência ao ataque sobre os silicatos hidrata- 
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tricálcico hidratado, mais facilmente do que-com os os 
silicatos: 


Ca (0H), + Na;S0,.10H, O = Ca S0,.2 H50 + 


+ Na 0H + 8 H,0 


2 (3 Ca0. Al,0, 12H, 0) + 3 (Na, SO, 10H, 0) — 


=3 Ca0. Al, 0, .3 CaSO, .31H,0 +2 AI (0H), 


+ 6 Na 0H + Ag 


3 Ca0.SiO,. Aq + Na, S0,.10H,0 = CaSO, .2H,0) + 


+ SiO, Na, +Aq 


Contudo, estas possíveis actuações dos sulfatos al- 
calinos sobre os compostos integrantes do cimento hidra- 
tado são relativamente pouco intensas e muito pouco no- 
civas para o betão, a menos que o equilíbrio das reac- 
ções seja continuamente deslocado, da esquerda para a 
direita, por eliminação do silicato de sódio, ou do alumi- 
nato de sódio, dos produtos da reacção. 

São principalmente os sulfatos alcalino-terrosos, 
de cálcio e de magnésio, que possuem maior poder de 
agressividade para o betão. Especialmente o sulfato de 
magnésio, que é muito solúvel e ionizável, é muito 
corrosivo para os silicatos e aluminatos hidratados, os 
quais se decompõem, segundo um esquema idêntico 
à acção dos sulfatos alcalinos, mas muito mais in- 
tensamente: 


3 Ca0.2SiO, . Aq + MgSO, . 7H, O = CaSO, . 2H, 0 + 
+ Mg (0H), + SiO, . Aq 
2 (3 Ca0.Al,0 


4 - Aq) + 3 (MgSO,.7H,0) = 


=3 Ca0./A1,0,.3 CaS0,.31H,50 + 3 Mg (0H), + 


Al (0H); + Aq 


Todavia, neste caso, o sulfoaluminato de cálcio 
em presença duma solução de sulfato de magnésio 
torna-se instável, e o equilíbrio da reacção é continua- 
mente desfeito pela sua acção sobre o hidróxido de 


cálcio resultante da decomposição deste sulfoalumi- 
nato; 
3 CaO.A1,0 -3 CaS0, .31 H,0 = Ca(0H), + 


+ 3(CaSO, 2H,0) + A | (0H) 3 +, + MgSO, . 7H,0 
3 Mg (0H), + 3 CaS0,.2H,0 + aq 


decomposição esta que se explica pelo facto de o pH 
da solução saturada de Mg(0H), ser inferior ao pH re- 
querido para o equilibrio do sulfoaluminato. 

De igual modo a decomposição dos silicatos de 
cálcio é uma consequência do mais fraco pH duma solu- 


89 


id 


ção saturada de hidróxido de magnésio (pH = 10,5) em 
relação ao pH = 12,4 da solução saturada de hidróxido 
do cálcio que mantém a estabilidae alcalina do be- 
tão [1]. 


(Nota; no caso da acção da água salgada o mecanismo 
químico da actuação do sulfato de magnésio sobre os 
aluminatos do cimento hidratado não se efectua como 
anteriormente se indicou, porque em presença duma 
solução de cloreto de sódio (= 3%) o aluminato 
tricálcico se decompõe em bicálcico que permanece es- 
tável; ou então combina-se com o NaCl, sob a forma 
de cloroaluminato de cálcio, se o concentração de sal 
for superior à indicada, o que parece ser o caso do 
betão imerso na água do mar). 

— Acção dos oxidrilos e carboxilos: Os produ- 
tos orgânicos portadores de oxidrilos ou de carboxilos — 
tais como a glucose, sacarose, fenóis, ácidos gordos, 
etc. — actuam fâcilmente sobre o hidróxido de cálcio 
resultante da hidratação do cimento formando os corres- 
pondentes derivados orgânicos do cálcio (alcolatos, 
sacaratos, fenolatos, gliceratos, sais de cálcio orgânicos, 
etc.). Na sua forma mais geral estas reacções com os 
produtos orgânicos, hidroxi ou carboxilados, pode esque- 
matizar-se: 


R-—-O 
2 R- OH + Ca l0H),=p . q>Ca + 2 H,0 


R — CO.O 


Tanto os ensaios laboratoriais como a prática têm 
demonstrado que a maior parte dos produtos orgânicos, 
químicos ou alimentares, não tem uma acção muito 
corrosiva sobre o betão. No entanto há certos ácidos 
gordos de óleos vegetais, por exemplo, que em condi- 
ções favoráveis podem desagregar o betão normal, e 
tanto mais intensamente quanto mais oxigenados ou 
insaturados forem. Parece que certos grãos que germi- 
nam em terrenos húmidos contêm ácidos orgânicos 
que podem destruir os melhores betões, pelo que, 
possivelmente, afectarão também os próprios silicatos 
e/ou aluminatos hidratados do ligante hidrálico [2]. 

Ainda no que respeita aos mecanismos gerais de 
ataque químico do betão não devem deixar de referir-se 
certos outros processos de corrosão que podem conduzir 
a expansão, perda de resistência, ou desagregação das 
estruturas de betão; é o caso, por exemplo, do ataque e 
destruição do aço das armaduras por agentes químicos 
e biológicos que penetram no interior do betão; ou o 
ataque dos agregados de natureza calcária por acção 
dos ácidos; ou ainda as reacções entre certos agregados 
vitrosos e os alcalis do cimento, etc., mecanismos estes 
que são porém pouco significativos para o caso que mais 
nos ocupa agora. 


3 — O ATAQUE DO CIMENTO POR PRODUTOS DE 
USO INDUSTRIAL 


Como seria de esperar, a estabilidade quimica rela- 


tivamente fraca do betão de cimento portland em muitos 
meios agressivos tem sido a causa motora de muito 
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trabalho ce investigação na procura de cimentos espe- 
ciais quimicamente mais resistentes para certos 
fins. Por outro lado, dirigiu-se também a atenção para 
a selecção dos agregados, dosagem e aplicação do betão 
em obra, utilização de aditivos convenientes para me- 
lhorar as características de fluidez, endurecimento, capa- 
cidade e impermeabilidade do betão, etc. Pelo que res- 
peita aos cimentos portland especiais há muito que se 
estabeleceram normas quase internacionais que, como se 
sabe, regulam, qualitativa e quantitativamente, o teor 
dos compostos integrantes do clinquer para determinadas 
aplicações do cimento em meios agressivos específicos 
— cimentos portland de alta resistência inicial, para 
resistência aos sulfatos, de baixo calor de hidratação, 
cimentos aluminosos, cimentos portland pozolânicos, etc. 


No que respeita ao assunto que agora nos interessa, 
um dos primeiros estudos experimentais, na Europa, 
sobre a decomposição dos cimentos por produtos de 
uso industrial foi publicado por G. Baire, e nele se sa- 
lientou a vantagem de utilizar o cimento pozolânico, em 
lugar do cimento portland normal [3]. Neste estudo o 
autor conclui «que o cimento pozolânico («ciment à la 
gaise), mesmo em caso de decomposição, apresenta 
sempre uma resistência aos produtos químicos apreciã- 
velmente superior ao cimento portland normal». O tempo 
de conservação ao ar, entre a amassadura e a imersão 
no reagente, tem grande influência no ataque, natural- 
mente para igual concentração das soluções agressivas. 
Os produtos ensaiados foram agrupados em três catego- 
rias, de acordo com a sua acção quimica sobre os dois 
tipos de cimento nas condições do ensaio (cubos de 
argamassa plástica a 1:3 imersos nos líquidos até me- 
tade da altura): 


(a) Os que atacam apreciâvelmente as duas espé- 


cies de cimentos: óleo de linhaça — óleo animal de 
bovinos — cloreto de amónio — sulfato de amónio — sul- 
fato de ferro (ataque fraco) — cloreto de magnésio em 


solução concentrada; 

(b) Os que atacam o cimento portland mas não o 
cimento pozolânico: açúcar cloreto e hipoclorito 
de cálcio (cloreto comercial) — glicerina — nitrato de 
amónio — sulfato de cobre — nitrato de sódio — cloreto 
de cálcio — cloreto de magnésio — sulfato de magnésio 
(mesmo em solução fraca) — sabão — hipossulfito de 
sódio -— sulfito de sódio — bicromato de potássio; 

(tc) Os que não atacam apreciâavelmente nenhum 
dos cimentos nas condições do ensaio: óleos minerais 
— Óleo de rícino — petróleo, benzol, e hidrocarbonetos 
em geral — álcool — soda cáustica não muito concen- 
trada — ácidos orgânicos em solução diluída (oxálico, 
cítrico, tartárico, tânico) — glucose — alúmen de po- 
tássio — sulfatos de sódio e de zinco — carbonatos 
alcalinos — sulfidrato de amónio — água de cloro — 
alcatrão e «mazout» — bissulfito de sódio, 


Já se disse anteriormente que um dos mecanismos 
de ataque do betão pelos produtos com ele em contacto 
pode ser desencadeado pela acção do elemento nocivo 
sobre a cal do cimento, a qual se desloca facilmente por 


iões tais como SO,”, CI , OH —, COOH, etc. Como 
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se sabe, no cimento pozolânico, esta cal (hidróxido de 
cálcio) está em proporção bastante mais fraca e «bloquea- 
da» pela sua combinação com a sílica solúvel da pozo- 
lana activa. Logo, como é natural, o betão com cimento 
pozolânico é mais estável ou menos vuinerável a este 
tipo de corrosão química, 


4 — PROTECÇÃO DOS DEPÓSITOS DE BETÃO. 
REVESTIMENTOS 


Para se evitar, ou atenuar, a corrosão química do 
betão, vários meios têm sido propostos e são praticados, 
desde o tratamento superficial com indutos silicatados 
(silicato de sódio, silicato de etilo, fluossilicatos de 
sódio, magnésio e zinco, etc.) ao uso de películas del- 
gadas de tintas e vernizes mais quimicamente resisten- 
tes, até aos revestimentos com barreiras mais espes- 
sas de produtos cerâmicos e argamassas cimento resina, 
ou com folheados de materais pláticos e metálicos. Po- 
rém, antes de se pensar em escolher o processo de re- 
vestimento mais conveniente em dado caso prático, é 
importante ter a certeza de que o betão precisa real- 
mente de ser protegido. Há, como vimos, muitas subs- 
tâncias que são inertes para com o betão, ou não o 
atacam apreciâvelmente, entre as quais se contam, como 
vimos, a malor parte dos carbonetos e nitratos, alguns 
cloretos e fluoretos e os silicatos; as soluções alcalinas 
fracas; os produtos derivados do petróleo desde que 
isentos de aditivos acídicos; e alguns destilados do 
alcatrão de carvão. Basta uma fabricação cuidada do 
betão para conseguir boas condições de armazenagem — 
qualidade e dosagem do cimento, granulometria conve- 
niente e boa estabilidade dos agregados, uso de certos 
aditivos melhorantes do betão, compactação e tratamento 
adequado após a colocação em obra, etc. Há outros casos 
em que a agressividade do produtos para com o betão 
é relativamente diminuta e será suficiente uma simples 
aplicação de fluossilicatos para penetrarem através do 
betão e se combinarem com alguma cal livre, formando 
fluossilicatos de cálcio insolúveis. Portanto, a primeira 
condição para escolher um revestimento do betão é que 
seja adequado à resistência química necessária na expo- 
sição requerida. Aliás, além da natureza química do pro- 
duto a armazenar outros factores podem determinar a 
agressividade para com o betão: a temperatura, a pres- 
são, as variações ciclícas de humidade e molhagem, 
assim como, evidentemente, a já mencionada qualidade 
e colocação do betão em obra. 

O grau de protecção dos depósitos, ou de quaisquer 
outras estruturas em betão, aconselhável numa dada cir- 
cunstância, dependerá assim principalmente do seguinte: 


— Natureza do produto ou agente químico. 

— Temperatura, concentração e pressão. 

— Volume do produto por unidade de superfície 
do betão. 

— Exposição contínua ou intermitente. 

— Características do betão fabricado e posto em 
obra. 


Os dados bibliográficos mais recentes sobre a 
protecção do betão contra ataque químico por meio de 
revestimentos orgânicos, e de outros materiais resisten- 
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tes à corrosão, foram publicados pelo «ACI Commitee 
515» sob a forma de manual, em Dezembro de 1966, 
pelo «Jornal of the American Concrete Institut» [4]. 
Nele se referem os diferentes métodos para a protecção 
anticorrosiva do betão e as espécies de revestimentos 
mais adequados à sua defesa contra a acção de cada 
produto ou reagente químico. Salienta-se, especialmen- 
te, a tabela final deste Manual-Guia na qual se indica, 
para cada produto ou reagente químico, o efeito que 
provoca sobre o betão normal e os respectivos tratamen- 
tos superficiais ou barreiras espessas mais aconselhá- 
veis em cada caso. 


Os diferentes materiais e métodos para proteger 
o betão são indicados resumidamente, como segue. 


4.1 — Tratamentos superficiais com películas inorgã- 
nicas silicatadas 


Esta técnica clássica de tratamento superficial do 
betão com indutos ou soluções de silicato de sódio ou 
de fluossilicatos de sódio, zinco, magnésio e alumínio, 
é um tanto eficaz, e pode prestar bons serviços em 
certos casos de menos responsabilidade em que não haja 
agressividade acentuada do produto para com o betão 
nem para com a própria película silicatada. O tratamen- 
to tem como principal efeito aumentar a dureza e im- 
permeabilidade das paredes do betão devido à reacção 
entre os compostos de cálcio solúveis do cimento hi- 
dratado (principalmente do hidróxido de cálcio) e a 
sua transformação em silicatos ou fluossilicatos de cál- 
cio insolúveis. Para melhorar o efeito destes produtos 
sobre o betão convém, portanto, que penetrem nele o 
mais profundamente possível para reagirem com a cal 
através duma camada espessa de betão, e não forma- 
rem simplesmente uma película superficial que venha 
a destacar-se ou a fender facilmente e ponha o betão 
a descoberto [5]. Por isso, a primeira camada deve dar- 
-se muito diluída para não fechar completamente os poros 
e capilares do betão, e as camadas subsequentes progres- 
sivamente mais concentradas. Estas soluções de silicato 
de sódio, ou de fluossilicatos de magnésio e zinco, não 
devem ser aplicadas ao betão muito fresco. Recomen- 
da-se cerca de 28 dias de idade, e a aplicação de 2 a 
3 camadas (0,120 kg por litro, de fluossilicato, para a 
1.º demão e cerca de 0,25 kg por litro para a camada 
final). 


O chamado «processo Ocrate» consiste no trata- 
mento do betão pré-moldado, com tetrafluoreto de silício 
gasoso (F, Si), numa camãra de vácuo na qual é introdu- 
zido o gás sob pressão. Nestas condições parece que 
o F, Si penetra mais profundamente no betão e o torna 
química e bacteriolôgicamente resistente, Aliás, os fluos- 
silicatos alcalinos e alcalino-terrosos são também bons 
insecticidas, o que pode ser benéfico em certos casos 
tal como na armazenagem de cereais e produtos análo- 
gos. Mas nunca se deve usar, nestes casos de contacto 
com produtos alimentares, o fluossilicato de chumbo 
porque é tóxico; o mais indicado, para este fim, será o 
fluossilicato de alumínio ou de sódio. Como já se disse, 
reduzem apreciâvelmente a penetração e absorção de 
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água através do betão mas, no entanto, são pouco efi- 
cientes contra a penetração de líquidos leves tais 
como as gasolinas e solventes análogos, O citado ma- 
nual da ACI aconselha a sua utilização no tratamento 
superficial do betão em contacto com os seguintes ma- 
terais e produtos químicos: 

— Ácido acético (10%) — Sulfato de alumínio — 
Cloreto de amónio — Gases de escape de automóveis 
(N, O, CO, NOs, SO5, 05, H50, vapor, hidrocarbonetos 
não queimados, etc.) — Hidróxido de bário — Gordura 
animal de bovinos — Benzeno — Hidróxido de cálcio — 
Nitrato de cálcio — Carbazol — Anidrido carbónico e 
óxido de carbono — Óleo de rícino — Criseno — Óleos 
de alcatrão de carvão — Óleo de algodão — Creosote 
— Cresol — Cumeno — Destilados de lavandaria — Clo- 
retos e sulfatos de ferro — Óleos de peixe — Sumos de 
frutos — Glucoses — Glicerina — Ácido húmico — 
Ácido sulfídrico- lodo — Ácido láctico — Soluções 
de refinação do chumbo — Óleo de lenhite — Óleo de 
linhaça — Óleos de lubrificação — Cloreto e sulfato de 
magnésio (não muito concentrados) — Sulfato de man- 
ganês — Estrumes — Mosto em fermentação — Cloreto 
mercuroso e mercúrico — Refugo de água de minas — 
Óleo mineral — Terebentina e «white spirity — Melaços 

Óleo de mostarda — Sulfato de níquel — Ácido oleico 
Azeite — Parafina — Fenol e fenantreno — Ácido fosfó- 
rico — Óleo de papoila — Sulfato de alumínio e potássio 
Dicromato de potássio — Sulfato e persulfato de 
potássio — Óleo de colza — Água do mar — Forragens 
Sais de sódio (excepto o hidróxido, cianeto, hipo- 
clorito, nitrato e monofosfato) — Óleo de soja — Açúcar 
— Resíduos do sulfito — Óleo de sebo — Ácido tânico e 
resíduos tanantes — Tabaco — Óleo de madeira da china 
— Óleo de noz — Óleo de baleia — Xileno — Cloreto 
e sulfato de zinco. 

Estes tratamentos superficiais do betão com produ- 
tos fluossilicatados previnem bastante o ataque pelos 
produtos citados; aliás, como seria de esperar, todos 
eles de agressividade pouco acentuada para com o ci- 
mento hidratado do betão. Por outro lado, podem ser 
usados com vantagem na preparação prévia da super- 
fície do betão para aplicação de certas tintas conven- 
cionais. Mas esta técnica de pintura está já um tanto 
ultrapassada pois, como veremos a seguir, existem hoje 
no mercado produtos de pintura, ou similares, que se 
podem aplicar directamente sobre o betão, esteja ele 
como estiver, sendo apenas aconselhável alisar a superfi- 
cie a revestir. 


4.2 Revestimentos orgânicos (tintas e vernizes) 


As tintas e vernizes usadas sobre o betão para evitar 
o seu ataque quimico serão, evidentemente, as de carac- 
teríticas antialcalinas, e capazes de resistir à acção 
quimica do agente agressivo com que vão ficar em con- 
tacto em cada caso específico. Realmente, dispõe-se de 
várias soluções possíveis consoante a importância e a 
qualidade do problema a resolver — água potável ou 
não, produtos alimentares sólidos ou líquidos, ácidos 
fortes ou fracos, solventes hidrocarbonados, produtos 
químicos diversos, etc. Estas soluções podem ir desde 
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os clássicos revestimentos termoplásticos com tintas 
antialcalinas e antiácidas do tipo betuminoso ou basea- 
das em borracha clorada, até às pinturas aquosas do 
tipo vinílico e acrílico, ou aos revestimentos especiais 
do género uretano, poliester e epoxílico de endureci- 
mento por catalisação a frio (geralmente a dois com- 
ponentes separados). Na protecção de depósitos de betão 
para armazenar solventes orgânicos, por exemplo, um dos 
sistemas de pintura que se revelou mais eficaz foi basea- 
do em copolímeros de cloreto de vinilideno e nitrilo 
acrílico, aplicado sobre uma demão de primário de 
borracha polissulfureto, A dispersão aquosa dos dois 
produtos dá excelentes resultados e elimina o possível 
transtorno de certos solventes voláveis das tintas quan- 
do se aplicam em locais de pouco arejamento; mas, em 
contrapartida, esta dispersão aquosa parece ser relati 
vamente pouco estável (cerca de 2 meses) [5]. Também 
o problema da protecção, por meio de pintura, dos depó 
titos de betão para armazenagem de vinhos e produtos 
similares é particularmente delicado devido a ser difi- 
cil garantir a absoluta inoquidade dos produtos, que 
não podem comunicar cheiro ou sabor de qualquer es- 
pécie ao conteúdo. Quase todas as pinturas usuais de- 
vem ser postas de parte, justamente porque transmitem 
mau gosto ao vinho [6]. O que mais se recomenda, neste 
género de protecção, parece ser a pintura com emulsões 
aquosas de copolímeros vinílicos (que devem deixar-se 
endurecer completamente antes de encher os reservató- 
rios) ou então as pinturas mais recentes do tipo epoxi- 
lico. Aliás estes revestimento orgânicos baseados em 
resinas de epoxilina catalisadas a frio, modificadas ou 
não com pez de alcatrão ou com polissulfuretos, são dos 
mais quimicamente resistentes, e por isso dos mais 
adequados à generalidade dos agentes agressivos para 
o betão, incluindo solventes orgânicos. Embora estas 
tintas já possam, como se sabe ser preparadas e forne- 
cidas em um só componente pronto a usar (pré-políme- 
ros bloqueados), são ainda mais vulgarmente comerciali- 
zadas sob a forma de dois componentes separados para 
se misturarem na altura da aplicação. Como se sabe, 
há diversos compostos que se utilizam como catalisa- 
dores do endurecimento destas resinas, a maior parte 
dos quais são compostos aminados (poliaminas e polia- 
midas) que podem apresentar certa toxicidade. Sempre 
que se trate de revestir depósitos de betão para arma- 
zenagem de produtos de natureza alimentar é indis- 
pensável ponderar este problema, e, de qualquer modo, 
proceder sempre ao tratamento da película com vapor 
de água; primeiro por uma questão de prudência contra 
qualquer eventual toxicidade transmissível pelo catalisa- 
dor aminado, e em segundo lugar porque o endureci- 
mento é beneficiado, adquirindo a película melhores 
qualidades de resistência em face do betão e do pro- 
duto a armazenar. Às indicações dos fabricantes ou for- 
necedores destes revestimentos baseados em resinas de 
epoxilina (tintas e vernizes «epoxy»), que devem ser 
seguidas criteriosamente, proclamam a sua assinalável 
resistência química tanto à água como a soluções sali- 
nas diversas e a ácidos e bases não muito concentrados. 
Igualmente resistem à acção de óleos e gorduras vege- 
tais e animais, à maioria dos solventes e carburantes, etc. 
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A sua aderência directa ao betão, e a outros materiais 
baseados em ligantes hidráulicos, é excepcionalmente 
boa, desde que se respeitem as indicações do fabri- 
cante no que se refere a condições de aplicação. Geral- 
mente a vida útil depois de se juntarem os dois com- 
ponentes «pot life» anda à volta de 15 a 20 horas. O 
tempo de secagem ao tacto (fora de poeiras) é cerca de 
1 hora e o intervalo para aplicar segunda demão deve ser 
de pelo menos 24 horas. O endurecimento completo faz- 
-se em 7 a 10 dias, consoante as condições ambientes, 
mas para serviço debaixo de água convém esperar até 
cerca de 15 dias. A película resiste bem a temperaturas 
de 80'C em calor seco, mas em calor húmido não aguen- 
ta muito tempo a esta temperatura, À aplicação do re- 
vestimento deve fazer-se fora da acção directa do sol 
e da chuva, de preferência durante as primeiras horas 
da manhã ou ao fim da tarde. Como se sabe a tinta 
contém solventes voláteis e inflamáveis, por isso deve 
trabalhar-se em locais bem arejados, ou então os operá- 
rios devem munir-se das respectivas máscaras de protec- 
cão. Para serviço em futuro contacto com água potável 
ou qualquer produto de natureza alimentar aconselha-se 
sempre medidas especiais para o endurecimento rápido 
e total eliminação de possivel toxicidade, tal como o 
tratamento da película seca pelo vapor de água, já antes 
mencionado, ou qualquer outro que seja eficaz, Convém 
também prolongar a entrada em serviço efectivo até 
cerca de 4 semanas após a aplicação da tinta. 


(Nota: as características indicadas dizem respeito à 
pintura «epoxy» a dois componentes separados, da 
marca «Colmaso! SIKA», mas são também as geral- 
mente indicadas para outras marcas de tintas do mesmo 
tipo, Convém sempre, no entanto, que o utilizador se in- 
forme junto dos fabricantes ou fornecedores, não só para 
estas tintas como para outras de tipo especial). 


Referem-se aqui, mais especialmente, as pinturas do 
betão com tintas baseadas em resinas de epoxilina por- 
que parece ser este actualmente um dos tipos de tintas 
mais indicados, e de utilização mais geral e adaptada 
ao caso da armazenagem de produtos químicos e alimen- 
tares que agora mais directamente nos ocupa. No en- 
tanto, deve acentuar-se que haverá muitos casos de 
menor responsabilidade em que outros sistemas ce pin- 


tura com as tintas usuais a solvente aquoso — tais 
como as emulsões betuminosas, viniílicas e acrílicas, e 
estireno — butadieno, ou com as tintas a solvents or- 


gânico baseadas no breu de alcatrão, em resinas fenó- 
licas, em borracha clorada, em polietileno, em silicones, 
ou noutros polímeros de tipo termoplástico, quer sim- 
ples quer modificados, serão perfeitamente adaptados e 
econômicamente mais favoráveis. É evidents que, como 
já se disse, tudo depende das circunstâncias específicas 
do caso prático a resolver e em função da agressividade 
respectiva, 

Como indicação sumária apontam-se a seguir al- 
gumas circunstâncias, embora relativamente poucas, em 
que não é viável, ou pelo menos não é segura, a pro- 
tecção do betão contra a corrosão por meio do revesti- 
mento com qualquer espécie de películas orgânicas 
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delgadas (Vid. «ACI Guide for the protection of concrete 
against chemical attack» [4]: 


— Acetona — Ácido crómico 50%, Ácido ni- 
trico, 40% — Ácido fluoriídrico, 75% — Ácido sulfúrico 
conc. — Adubos — Dinitrofenol — Sulfureto e tetraclo- 
reto de carbono — Metiletilcetona e metilisobutilcetona 
— Tricloroetileno — Óleos de alcatrão de carvão — Gases 


quentes — Extractos fermentados de frutos, vegetais 
e cereais — Sulfato de zinco. 
4.3 — Revestimentos cerâmicos e de vidro. Argamassas 


quimico-resistentes 


Os tijolos especiais e os azulejos cerâmicos ou de 
vidro são, de longa data, usados no revetimento do betão 
contra o ataque químico, quando seja necessário formar 
uma barreira bastante mais espessa do que as películas 
superficiais orgânicas e inorgânicas. Alguns tipos destes 
produtos cerâmicos são especialmente fabricados para 
o efeito, e a ASTM C-279 especifica dois destes tipos 
para usar em função da absorção e do choque térmico. 
Na actualidade, a utilização dos revestimentos cerâmicos 
do betão tem progredido acentuadamente devido ao 
desenvolvimento de argamassas para assentamento do 
material que possuem melhores caracteríticas de resis- 
tência ao ataque pelos agentes agressivos do que as nor- 
mais argamassas hidráulicas. Por outro lado, o problema 
do isolamento das juntas, muito importante na maior 
parte dos casos, tem também melhorado considerável- 
mente com o aparecimento no mercada dos novos elas- 
tôómeros e polímeros sintéticos para preparação dos 
tapa-juntas e a sua fácil aplicação com «seringa» [7, 11]. 

Os principais tipos de argamassas colantes hoje 
usadas nas alvenarias quimico-resistentes são os seguin- 
tes: (a) as argamassas de enxofre plastificado com 
polissulfuretos que são geralmente apresentadas na for- 
ma de lingotes e aplicadas por fusão a quente; (b) as 
argamassas de silicato, há muito usadas em barreiras 
protectoras quimico-resistentes, as quais consistem numa 
solução do silicato de sódio e «fillery misturado com um 
agente endurecedor ácido que precipita a sílica-gel 
quando se junta ao silicato (geralmente trata-se de um 
sílico-fluoreto); (Cc) as argamassas com resinas fenólicas 
modificadas e de furano, adicionadas de «fillers» diver- 
sos e endurecidos com catalisador ácido; (d) as arga- 
massas de resinas poliester insaturadas, adicionadas de 
estireno e catalisadas com um peróxido na altura da 
aplicação; (e) as recentes argamassas quimico-resisten- 
tes constituídas por resinas de epoxilina líquidas, mis- 
turadas com um «filler» e endurecidas com aminas alifá- 
ticas ou, mais frequentemente, com poliamidas. Estas 
argamassas de resinas «epoxy» têm propriedades físicas 
e mecânicas excelentes, boa resistência quimica aos 
ácidos e alcalis, e podem ser preparadas formulações 
especiais para resistirem aos solventes e ao calor [8]. 

Outros materiais são também usados como arga- 
massas quimicamente resistentes em casos especiais. 
Pode dizer-se que qualquer resina sintética na forma 
líquida, ou que possa ser dissolvida num solvente usual, 
pode servir de base a uma argamassa química. Como 
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se sabe, os materiais betuminosos, adicionados ou 
não de inertes, podem usar-se como argamassas resisten- 
tes aos ácidos e aos alcalis para serviço a baixa tem- 
peratura, aplicados como «mastic», ou por fusão a 
quente, ou ainda como folhas que podem conter ou não 
um elastómetro. Também, em muitas circunstâncias, a 
simples argamassa de cimento pode ser mais convenien- 
te, e mais barata, que as argamassas de enxofre, de sili- 
cato, ou de resinas sintéticas, tal como no caso de ser- 
viço em condições alcalinas moderadas que não afectam 
o ligante hidráulico. 

A prática do assentamento dos tijolos e azulejos 
para revestimentos cerâmicos sobre o betão com arga- 
massas químico-resistentes é recomendada nas normas 
da American Society of Testing Materials, especialmente 
as seguintes: 


ASTM C836 -— Use of sulfur-type chemical-resis- 
tant mortars. 

ASTM C397 — Use of chemically setting silicate- 

type chemical-resistant mortars. 

Use of hidraulic cement mortars in 

chemical resistante masonry. 

ASTM C399 — Use of resin-type chemical-resis- 
tant mortars, 


ASTM C398 


4.4 — Revestimento do betão com materiais folheados 


O revestimento das paredes de betão com materiais 
folheados, quimicamente resistentes, adquiriu maior 
versatilidade pelo uso de novos tipos de colas sintéticas 
capazes de colar o betão a não importa que tipo de 
material — borracha, plástico, madeira, metais, e ao pró- 
prio betão [9,10]. Contudo, a cola e a folha devem 
ser convenientemente escolhidas no que respeita aos 
líquidos a que possam ser expotas conjuntamente. É bas- 
tante importante evitar qualquer combinação que possa 
conduzir à difusão do solvente através da folha e à con- 
sequente solvatação da película de cola, 

Os materais usados em forma de folhas para pro- 
teger o betão contra ataque químico são as borrachas 
naturais e sintéticas, tais como o neopreno, cloreto de 
polivinilideno-nitrilo acrílico, borracha butilo, polissul- 
fureto, ou borrachas de polietileno clorossulfonado. No 
entanto o elastómero ainda mais usado parece ser a 
borracha natural, em folhas com a espessura desde 3 mm 
a 12,7 mm (1/8" à 1/27). A folha é normalmente vulca- 
nizada com vapor a alta pressão ou com água quente. 
As folhas de borracha de neopreno protegem bem o 
betão contra a maior parte dos ácidos, excepto o ácido 
clorídrico, e têm uma resistência considerável aos sol- 
ventes alifáticos. Para resistência a solventes é também 
muito útil a borracha de polissulfureto. A borracha de 
nitrilo e o polisobutileno encontram a sua principal 
aplicação na protecção do betão contra oxiácidos de 
concentrações moderadas; quanto à folha de borracha 
butilo possui excelente resistência ao ozono ou à difusão 
de gases. O polietileno clorossulfonado é largamente 
usado mas como revestimento de cobertura. 

No que respeita aos folheados de plástico, geral- 
mente reforçados com fibras de vidro, os mais usados 
para revestir o betão são as folhas de poliester, de clo- 


94 


reto de polivinilo e polivinilideno, os folheados de resinas 
epoxílicas, e os plásticos fluorcarbonados (tetraclorofluo- 
retileno e triclorofluoretileno). Estes últimos são resis- 
tentes a quase todos os produtos químicos, excepto me- 
tais alcalinos e fluor, aguentando mesmo a acção do 
ácido nítrico concentrado a temperaturas de cerca de 
205"C (400'F). Mas o seu uso como revestimento pro- 
tector do betão é um tanto limitado pela vulnerabilidade 
das colas requeridas para a sua colagem. 

Quanto a folheados metálicos o mais usado é a 
folha de chumbo resistente à corrosão (99,90% Pb; 
0,04 a 0,08% Cu; 0,002 a 0,20 Ag, menos que 0,005% 
Bi) que funde a cerca de 326 C e tem um peso especi- 
fico de 11,3. As folhas de chumbo do comércio têm a 
espessura de 0,4 mm até 6,3 mm (1/64 a 14”) e di- 
mensões superiores a 243 cm X 609 cm (8 X 20 ft) 
para minimizar o número de juntas sem que, no entanto, 
se tornem muito difíceis de manobrar, 

Na bibliografia consultada não se nos deparou qual- 
quer referência sobre a possível utilização de folheados 
de madeira e seus derivados no revestimento de depó- 
sitos ou outras estruturas de betão contra a agressivi- 
dade química. É certo que a madeira é um produto rela- 
tivamente poroso e que a sua utilização em folheados 
delgados não parece muito aceitável para servir de reves- 
timento anticorrosivo do betão, tanto mais que os fo- 
lheados de plástico são práticamente inertes e imper- 
meáveis. No entanto há certos aglomerados de madeira 
que se podem tornar também pràticamente impermeáveis 
— tal como o «Platex» nacional, por exemplo — e que 
são comercializados, como se sabe, em folhas de várias 
espessuras e dimensões. À primeira vista parece um 
material susceptível de substituir o revestimento do be- 
tão com folhas de plástico, nos casos práticos em que 
este método de protecção seja aconselhável, e com 
apreciável vantagem económica, À sua colagem ao betão 
será até mais fácil que a do plástico; e quanto à resis- 
tência a produtos químicos e alimentares, talvez mesmo 
a produtos em fermentação, parece um assunto que me- 
rece ser investigado experimentalmente. Tudo leva a crer 
que o referido aglomerado de madeira, ou outros do 
mesmo género, possa prestar bons serviços no revesti- 
mento de recipientes de betão para armazenar produtos 
alimentares em geral. O interesse será ainda maior caso 
a sua resistência se amplie aos produtos fermentados 
visto ser um dos casos em que se torna mais difícil 
encontrar processo económico e adequado de proteger 
as estruturas de betão por meio de revestimento. 


5 — RESUMO E CONCLUSÃO 


Devido ao uso cada vez maior do betão de cimento 
portland tanto no domínio das obras públicas como no 
das técnicas e operações industriais, incluindo a armeze- 
nagem de pordutos químicos e alimentares os mai diver- 
sos, estão continuamente surgindo novos problemas rela- 
tivos ao seu ataque químico, ao mesmo tempo que vão 
aparecendo no mercado novos produtos para a sua pro- 
tecção, cuja adaptabilidade a dadas condições de corro- 
são deve ser antecipadamente confirmada em cada caso 
prático específico. 
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Nos capítulos 2 e 3 do presente trabalho bibliográ- 
fico, fazem-se algumas considerações gerais sobre o 
ataque químico do cimento por produtos de uso indus- 
trial; e no capítulo 4 afloram-se os processsos mais viá- 
veis de proteger as superfícies de betão por meio de 
revestimentos anticorrosivos, a saber: 


(1) Tratamentos superficiais com produtos inor- 
gânicos siliciosos, especialmente silicatos e fluossilicatos 
de sódio, zinco, magnésio, e produtos análogos. 

(2) Revestimentos com películas de tintas e 
vernizes convencionais e especiais, particularmente as 
pinturas com produtos «epoxy» de endurecimento a frio. 

(3) Revestimento com uma barreira mais espessa 
de produtos cerâmicos ou de vidro — tijolos, azulejos, 
etc. — e com argamassas especiais quimico-resistentes. 

(4) Revestimento com folhas de plástico ou de 
metal (folhas de chumbo) coladas ao betão; ou talvez, 
em casos especiais, com folheados de madeira aglome- 
rada que tenham sido prêviamente tratados para se tor- 
narem práticamente impermeáveis. A este respeito parece 
aconselhável um estudo experimental sobre a resistên- 
cia química dos aglomerados de madeira aplicados como 
revestimento do betão, em especial na armazenagem de 
produtos químicos e alimentares, incluindo os produtos 
em fermentação. 


A decisão do método ou métodos mais aconse- 
lháveis deve ser ponderada, para cada caso, em face da 
natureza química do produto a armazenar ou em contacto 
com as paredes do betão; da temperatura, pressão e 
concentração desse produto; do volume por unidade 
de superfície em contacto; e da sua acção continuada 
ou intermitente, Certos produtos não atacam sensivel- 
mente o betão e a sua protecção torna-se desnecessária 
se ele for devidamente fabricado e posto em obra, Nou- 
tros casos a agressividade para o betão é relativamente 
pequena, de modo que bastará um tratamento da super- 
fície com um produto inorgânico endurecedor do cimen- 
to, geralmente um fluossilicato, seguido ou não duma 
pintura adequada às circunstâncias. Porém, no caso de 
agressividade acentuada dos produtos para com o betão 
é indispensável a sua protecção com uma película de 
pintura de tipo especial — mais geralmente epoxílica de 
endurecimento por catalisação a frio e tratada ou não 
pelo vapor, Ou então com uma barreira, mais espessa e 
impermeável, de produtos cerâmicos assentes em arga- 
massas resinosas químico-resistentes, ou com folhea- 
dos metálicos e de plástico convenientemente colados 
ao betão. Em qualquer dos casos é indispensável uma 
conveniente vedação com tapa-juntas resistentes ao meio 


agressivo. Em circunstâncias especiais — como seja, 
por exemplo, a armazenagem de produtos químicos ali- 
mentares — é evidente que tem de se entrar em conta 


com a possível alteração que esses produtos possam 
vir a sofrer no seu aspecto, gosto ou cor; além do pro- 
blema, essencial, da possível formação de toxicidade 
em face do betão e da película ou barreira protectora. 

Os dois pontos fracos do betão no ataque por agen- 
tes agressivos são o hidróxido de cálcio libertado na 
hidratação do cimento — e que torna o meio francamente 
alcalino, no qual os silicatos e aluminatos de cálcio se 
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estabilizam por cristalização — e, por sua vez também 
a fraca resistência química destes componentes hidra- 
tados e cristalinos que constituem o betão já endurecido. 
Portanto, independentemente da protecção que haja 
a fazer aos depósitos e outras estruturas de betão, por 
algum dos métodos de revestimento antes referidos, 
qualquer processo ou meio técnico no fabrico e na colo- 
cação do betão em obra que tenda a impedir uma ou 
outra destas causas da sua fragilidade química será 
benefício certo para a estabilidade da estrutura, em qual- 
quer circunstância ou meio ambiente em que tenha de 
prestar serviço. Especialmente quando se trate de con- 
tacto com líquidos leves que penetrem fácilmente através 
da superfície do betão, a agressividade será sempre mais 
reduzida, para iguais circunstâncias, se o betão for: 
(a) bem doseado e constituído por agregados quimica- 
mente estáveis com a granulometria adequada a uma boa 
compacidade; (b) não se expandir ou contrair, além do 
normal, para que não se formem fissuras ou fendas, du- 
rante o endurecimento, por onde possa penetrar mais fã- 
cilmente o líquido corrosivo; (c) for convenientemente 
posto em obra e devidamente cuidado no período a seguir 
à sua aplicação; (d) e, finalmente, lançando mão dos adi- 
tivos usuais aconselháveis para regularizar a presa e en- 
durecimento do cimento, e que protejam, ou fixem, o 
hidróxido de cálcio tanto no interior como à superfície 
do betão. 

Para definir e estabelecer o processo mais adequado 
de revestir o betão de cimento portland em face de cada 
produto industrial ou de natureza alimentar, aconselha-se, 
como ponto de partida, a consulta do livro «Guide for 
the Protection of Concrete Against Chemical Attack by 
Means of Coatings and Other Corrosion — Resistant 
Materials», publicado pelo ACI Committee 515, com es- 
pecial referência à sua tabela final e ao capítulo de pre- 
cauções e aplicação. 
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FORMAS MATRICIAIS NO CÁLCULO DE AFINORES (11) 


RESUMO 


Adaptam-se a algumas operações sobre afinores, 
as formas matriciais. 


1 — INTRODUÇÃO 


Neste artigo adaptaremos a algumas operações de- 
finidas no âmbito do cálculo de afinores, as formas 
matriciais, dadas as vantagens que advêm das proprie- 
dades do cálculo matricial. 

Consideraremos os afinores definidos num espaço 


vectorial topológico =,. 


2 — PRODUTO DE AFINORES 


O produto de dois afinores é um afinor cujas variân- 
cias são as dos afinores factores, sendo as componentes 
do afinor produto obtidas pelo produto das componentes 
dos afinores factores, de tal modo que os índices to- 
mam os mesmos valores nas componentes dos afinores 
factores e na componente do afinor produto 

Sejam os afinores 


É r 


12 AP JA 19 


z vezes contravariante 
vezes covariante 
T+ 5 V) 


f É 
2é AM JAM 


- vezes contravariante 


ivezes covariante 
di 
Y -— = Va 


O seu produto será o afinor. 
8 oo=— 92 AP+HIJA]a+m 
VE ENE Eis E spo Has 4 
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SYNOPSIS 


Matricial forms are adapted to some operations 
with afinors. 


a v vezes contravariante 
5 + Ívezes covariante 
a Pers 0 = Vi + Va Ya 


O produto de afinores goza das seguintes pro- 
priedades: 

— associatividade 

— existência de um elemento neutro, único, que é 
o afinor de ordem zero (escalar), em que a sua única 
componente toma o valor 1. 

Em cálculo matricial, as componentes do afinor pro- 
duto obtêm-se multiplicando directamente as matrizes 
coluna das componentes dos afinores factores: 


Fo 83) GE | ana Ea 


invl, 1) inve, 1) inv3, 1) 


3 — TRANSVECÇÃO DE AFINORES 


Transvectar dois afinores em relação a dois indices, 
um de cada afinor, de variância contrária, é somar os 
produtos dos dois afinores, dando ao par de índices a 
transvectar todos os valores de 1 an. 

O afinor resultado desta operação tem menos duas 
variâncias que a soma das variâncias dos afinores 
factores. 

Sem os afinores: 


4º Aejáa 


z vezes variante. 
2 AMAIm 


à vezes variante. 
A transvecção destes dois afinores em relação ao 
indice u, darã o afinor: 
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6º Ae +14] qt+m 
x + 3— 2 vezes variante 
Em cálculo matricial, as componentes do afinor 


resultado da transvecção, obtêm-se pelo produto das 
matrizes das componentes dos afinores factores: 


em que a arrumação das componentes dos afinores 
factores é feita dando aos z-1 (5-1) restantes índices 
os números de 2-1 (5-1) algarismos da base n, e 
colocando-os nas colunas (linha), para cada valor do 
indice a contrair. 


3.1 — Subiir e baixar indices de um afinor 


Esta operação é um caso particular de transvecção, 
em que um dos afinores é o tensor fundamental na forma 
contravariante ou covariante. 


ge SP |STa Jum — 


e por isomeria obtemos «: Ae|A4]/a9 
-m 
Em cálculo matricial, as componentes do afinor 


resultado da transvecção, obtêm-se pelo produto das 
matrizes: 


(n2—1 .n) 


(n,n) 


3.1.1 — Formas contravariante e covariante de um vector 


Como caso particular, teremos a correspondência 
entre as formas contravariante e covariante de um vector: 


Wk 9km Wm 


wi who wn | | W 


| Jim 


9nm 


(1,n) (n,n) (1,n) 


(1,n) (n,n) 


32 — Produto interno de vectores 


O produto interno de dois vectores é definido, 
em cálculo de afinores, pela dupla transvecção: 


U.v é = O mk * um = vk 


Em cálculo matricial, esta operação efectua-se do 
seguinte modo: 


PRP CE E v! | 8 
' 
ê 
(1,n) (n,n) (n,1) (1,1) 
3.3 — Norma de um vector 


À norma de um vector é definida, em cálculo 
de afinores, pela dupla transvecção: 


Jull = gn; - um. u 
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Em cálculo matricial, teremos, idênticamente: 
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APLICAÇÕES DO MÉTODO DE MONTE-CARLO AO 
ESTUDO DE FLUIDOS DENSOS 


RESUMO 


Apresentam-se em revisão algumas aplicações do 
método de Monte Carlo ao estudo de fluidos densos, 
abordando-se alguns dos modelos utilizados na simula- 
ção das propriedades dos fluidos densos, como os discos 
rigidos, esferas rigidas e Lennard - Jones; resume-se 
ainda a aplicação à simulação de cadeias poliméricas 
lineares. 


1 — INTRODUÇÃO 


O método de Monte - Carlo pode definir-se como 
o ramo das matemáticas que trata de experiências com 
números aleatórios. 

Embora algumas técnicas de Monte - Carlo muito 
simples fossem conhecidos desde o século passado, 
pelo menos, a aplicação sistemática do método, sobre- 
tudo a problemas de física nuclear e de investigação 
operacional, só começou a ser feita a partir do fim da 2.º 
Guerra Mundial(). 

O interesse deste método no estudo de fluidos 
densos resulta da possibilidade de resolução de integrais 
muito complexos. Na realidade, aplicando a Termodinã- 
mica Estatística, é possível, em princípio, calcular as 
propriedades físicas macroscópicas de um sistema qual- 
quer a partir das propriedades moleculares, em especial 
do potencial intermolecular. 

Para alguns casos, como gases diluídos e alguns 
sólidos cristalinos, este cálculo simplifica-se porque as 
particulas que os constituem podem ser consideradas 
como independentes, No entanto, para fluidos densos 
esta simplificação não é válida, tornando impossivel 
a resolução algébrica do problema. 

Para contornar esta dificuldade recorre-se a méto- 
dos computacionais que se podem classificar em duas 
classes: Dinâmica Molecular e Monte Carlo. 

O método da Dinâmica Molecular simula um siste- 
ma de N partículas e procede à integração das equações 
de Newton para o seu movimento simultâneo. 

No método de Monte Carlo estuda-se um sistema, 
constituído por N partículas, fazendo-o percorrer um per- 
curso aleatório no espaço configuracional, conhecido por 
cadeia de Markov, conforme se descreve adiante, 
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SYNOPSIS 


À review of some applications of the Monte Carlo 
methods to the study of dense fluids is presented. Some 
of the models used in the simulation of the dense fluid 
properties are discussed, such as the hard rod, the rigid 
sphere and the Lennard-Jones models. The application 
to the simulation of the chain molecules is presented 
as well, 


Em qualquer dos casos, o cálculo dos valores das 
grandezas física do sistema é feito com base no cálculo 
da energia potencial configuracional do sistema. É por- 
tanto necessário arbitrar um determinado potencial inter- 
molecular. 

Por outro lado, os dados fornecidos por estes mé- 
todos são, pela sua natureza, mais semelhantes a dados 
experimentais do que a resultados de cálculos teóricos, 
dado que resultam da simulação de um determinado sis- 
tema. Podemos portanto, considerá-los como verdadeiras 
experiências em computador sobre um sistema com um 
determinado potencial intermolecular. Estas experiências 
podem ser utilizadas de dois pontos de vista diferentes: 
ou se comparam os resultados obtidos com resultados 
experimentais obtidos para um determinado tipo de molé- 
culas, de modo a determinar a natureza do verdadeiro 
potencial intermolecular, ou então comparam-se com Os 
resultados obtidos teóricamente para um potencial inter- 
molecular, de modo a melhorar a teoria. 


2 — DESCRIÇÃO DO MÉTODO 


O método de Monte Carlo é directamente aplicável 
ao cálculo de qualquer valor expectável. 


fa fb) Px) dx (1) 
< í tg E né sa a E 
/ à P (x) dx 


em que P(x) é uma densidade de probabilidade não 
normalizada num espaço L-dimensional e f(x) uma 
função do vector x desse espaço, estendendo-se os inte- 


E 


grais à região !. 
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Para aplicarmos a expressão (1) a um sistema 
físico, precisamos de conhecer x e P(x). 

A Termodinâmica Estatística conduz-nos, dado um 
sistema constituído por N partículas, ocupando um vo- 
lume VW, à temperatura T, aos seguintes resultados: 


N (2) 


P(x) = exp — U (fr, ra, pa cs tg RT (3) 
em que U é a energia potencial configuracional do siste- 
ma, k a constante de Boltzmann, r, o vector de posição 
da partícula i e dv, o elemento de volume a ela asso- 
ciado. 

Para calcular o valor de <« f - dado por (1), utiliza- 
-se, no método de Monte Carlo, um percurso aleatório 
no espaço das configurações (espaço a 3N dimensões, 
das 3N coordenadas espaciais das N partículas que 
constituem o sistema). Parte-se de uma configuração ini- 
cial arbitrada e vão-se gerando sucessivamente novas 
configurações, de tal modo que configurações com a 
mesma energia potencial apareçam com uma frequência 
que tenda para o valor de P(x), dado por (3), quando o 
número de configuração geradas tende para infinito. Um 
processo deste tipo assegura-nos que a média de uma 
determinada propriedade avaliada sobre todas as n 
configuração geradas 


] n (4) 


tende para <« f - quando n -- O, 

A sucessão de pontos do espaço das configurações 
obtidas por este processo designa-se por cadeia de 
Markov (1.2), Para obtermos uma cadeia de Markov pre- 
cisamos de definir, além de uma probabilidade de ocupa- 
ção (ponto no espaço das configurações), a probabili- 
dade de transição entre um estado i e um estado j, P;; 

A Termodinâmica Estatística dá-nos a primeira 
probabilidade e o problema que se nos depara consiste 
portanto na definição duma matriz de probabilidade p ii 
que satisfaça à condição. 


lim f=<tf> (5) 


Para definirmos esta matriz de probabilidade, temos 
que definir as condições em que o problema se põe. 
Dada então uma configuração x, à qual corresponde 
um parâmetro t, dá-se origem a uma nova configuração 
x', à qual corresponde o parâmetro t + 1. 

Esta nova configuração está sujeita à condição de 
pertencer à vizinhança da configuração anterior x, defi- 
nindo-se vizinhança como o conjunto de configurações 
para as quais uma das N partículas, i, é deslocada da 


sua posição não variando cada uma das coordenadas 
de espaço mais que uma distância d arbitrada, e perma- 
necendo todas as outras na posição anterior. 

Os computadores digitais costumam representar as 
coordenadas de posição x através dum dado número 
de algarismos binários ou decimais, de tal modo que 
o espaço configuracional não é representado como um 
espaço contínuo, mas discreto, embora a densidade de 
pontos seja muito grande. 

Portanto à vizinhança da configuração y, « (y), 
pertence um número Y finito de configurações. 

Conhecido o potencial intermolecular, pois este é 
um dado do problema, pode calcular-se a energia po- 
tencial configuracional correspondente a cada configu- 
ração e daí obter o factor de Boltzmann dado por (3). 

À configuração x' é gerada com base nos números 
aleatórios: ,,:,,:4-4 -s, pertencentes ao intervalo [0,1]*. 

As partículas estão numeradas de 1 a N. Escolhe-se 
a partícula i a movimentar a partir de *,: 


= 1(N 4) (6) 


sendo | (x) a característica de x. As novas coordenadas 
desta partícula, (x', y, z), obtêm-se a partir das coorde- 
nadas iniciais, (x,y,z) por meio das equações (7), 
(8) e (9). 


x =x+d(1 -28) (7) 
V=y+dt -26) (8) 
*t=2+d1—-2Zhk (9) 


em que d define a vizinhança da configuração x, tal como 
foi descrito acima, 

Calcula-se a energia potencial da configuração x' 
e compara-se o valor de P(x')/P(x) com o número <.. 

Se “.<. P(x')/P(x), o estado x (t+ 1) é o novo 
estado x”. Se :c > P(x')/P(x), o estado x (t+ 1) é o 
estado inicial x, sendo t o parâmetro inteiro associado 
a cada configuração da cadeia de Markov. 

Percorre-se assim um percurso no espaço das con- 
figurações, aleatório, porque os números “, “5, Ey “4 
e *, são números aleatórios e uniformemente distribuidos 
no intervalo [0,1], 

A matriz de probabilidade de transição é 


=0 se jercti) 


j:rti)eP(j) > Pli) (10) 


"+ se 
13 


mem se | r (ij) e P(j) < P(i) 


Temos portanto como dados do problema uma configura- 
ção inicial do sistema de N partículas, ocupando um 
volume V, a uma temperatura T, e um determinado po- 


* NOTA: Estes números são na realidade pseudo-aleatórios, pois são gerados pelo computador através duma relação 
funcional com os números anteriores; procura-se sempre que esta função seja suficientemente complicada para que, para o 
efeito que se pretende, os números assim originados possam ser considerados aleatórios. 


102 


TÉCNICA 417 


tencial de interacção entre as particulas. Através do 
processo esboçado acima consegue-se determinar o valor 
expectável < [ - duma propriedade física cuja dependên- 
cia de V, T, N é dada por f. 

Para sistemas físicos comuns, N é pelo menos ca 
ordem de 1023, No entanto, para um número de parti- 
culas desta ordem de grandeza o tempo de computação 
e a capacidade de memória necessários, mesmo com 
os calculadores electrónicos modernos, tornam perfeita- 
mente impraticável a execução de um programa. 

Na realidade um sistema com N 1000 já seria 
muito difícil de tratar. Na tabela | apresentam-se alguns 
dados relativos ao tempo de computação, indicados por 
Me DONALD e SINGER 8. 


TABELA | 


Número de configu- 


ad coa ia rações por hora 
IBM 704 108 6 500 
ICT ATLAS 108 260 000 
UNIVAC 1108 864 400 000 


O que se pretende quando se utiliza o método de 
Monte Carlo é que os resultados reproduzam os valores 
de propriedades macroscópicas de sistemas fisicos e o 
facto de não podermos simular sistemas com um número 
de partículas da ordem do de sistemas macroscópicos 
reais é uma limitação que convém diminuir ao máximo. 
O principal erro introduzido pelo pequeno número de 
partículas é o da multiplicação dos efeitos de superfície, 
isto é, o número de partículas que não estão completa- 
mente rodeadas, até uma distância considerável compa- 
rada com o alcance do potencial intermolecular, por ou- 
tras partículas pertencentes ao sistema passa a ser pro- 
porcionalmente muito maior que nos sistemas macroscó- 
pico. Para evitar este efeito, METROPOLIS et al(4) 
imaginaram o artifício de considerar a célula de volume 
V contendo N partículas a uma temperatura T, rodeada 
por todos os lados por outras células idênticas, com o 
mesmo número de partículas, reproduzindo integralmente 
a configuração na célula fundamental.O sistema total 
passa a ser constituido pela célula fundamental e pelas 
«células imagem», em número infinito. Quando a par- 
tícula i da célula fundamental é movimentada pelo com- 
putador, as partículas |', à”, . . ., das «células imagem» 
têm um movimento exactamente idêntico. Para deter- 
minar as coordenadas destas basta portanto efectuar 
um deslocamento da origem do referencial. 

Os efeitos de superfícies são assim eliminados. 

A fig. 1 representa um sistema do tipo descrito 
acima, a duas dimensões. 

Com base num conjunto de células deste tipo, o 
cálculo das propriedades configuracionais macroscópicas 
do sistema é feito a partir das expressões obtidas pela 
Termodinâmica Estatística para essas propriedades, em 
função do potencial intermolecular. 

A energia potencial é calculada partindo da hipótese 
da aditividade, isto é, a energia potencial configuracional 
do sistema é tomada como a soma das energias de inte- 
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Fig. 1 — Método de Monte Carlo a duas dimensões. Configura- 
ção aleatória na célula central, rodeada por células-imagem 


racção entre todos os pares de particulas. Por outro lado, 
como a dimensão do sistema se tornou infinita, devido 
ao artifício das células imagem, o cálculo da energia 
potencial de cada uma das partículas na célula funda- 
mental torna-se impossível para potenciais que se anulem 
só para uma distância infinita. Para evitar esta dificuldade 
suplementar, a energia potencial calcula-se como a soma 
de duas parcelas, sendo a primeira relativa à interação 
entre cada uma das partículas da célula fundamental e 
as que se encontrem a uma distância menor que L, um 
comprimento adequado, e a segunda calculada por es- 
timativa. 

A uma aproximação deste tipo chama-se truncagem 
de potencial. 

Uma vez calculado o valor médio da energia con- 
figuracional, “U- a energia interna é dada por: 


3 (11) 
E= ss NKT+<US 
2 
A equação de estado é dada pelo teorema do virial 
pV <W > (12) 
NkKTOo 3NkT 
em que W é a função do virial 
I dulr;;) (13) 
W = 8 fa 
2 “dr 


sendo r;; e a distância entre as partículas i e je u( r;;) 
o potencial de interacção. O símbolo no somatório 
indica que jHi. 

Outra função calculável pelo método de Monte 
Carlo é o calor específico a volume constante 


(14) 

Roe RA aê U>? 
E —-—— Nk == + P — <US* 
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A função de distribuição radial, g(r), função de 
grande importância no estudo das propriedades configu- 
racionais de líquidos também pode ser calculado por 
Monte Carlo, mas a dedução da expressão que permite 
o seu cálculo reveste-se de tais dificuldades matemáticas 
que preferimos não a apresentar aqui. 


3 — LIMITAÇÕES DO MÉTODO 


A limitação principal da aplicação do Método de 
Monte Carlo ao estudo termodinâmico dos fluidos densos 
resulta de não ser possível calcular propriedades de 
transporte. 

No entanto, quando aplicado ao cálculo de valores 
de propriedade de equilíbrio, o método de Monte Carlo 
é exacto, isto é, não está sujeito a nenhum erro sistemá- 
tico, desde que seja possivel evitar os chamados proble- 
mas ergódicos, de que falaremos adiante, 

Evidentemente, existe sempre um determinado erro 
estatístico, Como já vimos na secção 2, f> constitui o 


limite de f só quando o número de configurações tende 
para infinito. Como o número de configurações é forçosa- 
samente finito, o valor calculado de <«f- está sempre 
afectado de um erro estatístico. Por outro lado, a necessi- 
dade de diminuir este erro, conjugada com as limitações 
da capacidade da memória e rapidez de cálculo do com- 
putador, implica, como já vimos, que o número N de 
partículas do sistema seja pequeno. O desvio em relação 
aos valores das propriedades de sistemas macroscópicos 
introduzidos por este facto pode ser estudado fazendo 
variar N e, segundo WOOD([2), as tentativas efectuadas 
parecem indicar que ele é proporcional a N”1, 

Quanto aos problemas ergódicos, eles consistem 

bâsicamente na existência de barreiras de potencial infi- 
nito ou muito elevado entre algumas regiões do espaço 
das configurações. Arbitrada então uma configuração ini- 
cial, a cadeia de Markov irá gerar novas configurações 
pertencentes sempre à mesma região do espaço das con- 
figurações, porque essas barreiras de potencial impedem 
a passagem para outra região. 
Suponhamos um sistema constituido por N esferas rigi- 
das numa configuração cúbica de faces centradas com- 
pacta. Se considerarmos uma das esferas numa posição 
fixa, temos no espaço das configurações (N—1)! pontos 
distintos, correspondentes às permutações das outras 
N-—1 esferas. Se arbitrarmos uma determinada confi- 
guração inicial, não é possível por meio de uma cadeia 
de Markov, atingir qualquer outra das configurações 
possíveis, porque para isso teriamos ce deslocar as esfe- 
ras de modo a penetrarem umas nas outras, O que corres- 
ponde a um potencial infinito. 

Este sistema constitui um exemplo extremo co que 
afirmamos acima, porque a sua densidade é a mais ele- 
vada possível. 

No entanto os problemas ergódicos aparecem sem- 
pre para densidades razoávelmente elevadas, quando a 
distância média entre o centros das partículas é inferior 
a duas vezes o diâmetro de cada uma delas. 

O facto de eles aparecerem vai afectar o processo 
de determinação dos valores médios das propriedades do 
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sistema, porque ao arbitrarmos uma configuração inicial, 
obrigamos o sistema a permanecer sempre na mesma 
região, e o valor médio não é calculado sobre todos 
os estados possíveis do sistema, mas só sobre os 
correspondentes à região em que o colocámos, o que 
limita considerávelmente o significado desse valor médio 
como valor duma propriedade física macroscópica. 


4 ALGUMAS APLICAÇÕES 
4.1 — Discos Rigidos 


A primeira aplicação do Método de Monte Carlo 
ao estudo dos fluidos densos deve-se a METROPOLIS et 
alt4) O trabalho destes autores consiste na determinação 
ca equação de estado, p = f (W,T,), de um sistema ce 
discos rígidos, que é um sistema a duas dimensões, 
constituído por partículas circulares, de diâmetro c e 
tomando o potencial intermolecular a forma 


o q G 
u tr) = (15) 
O, Ff v) 


em que r designa a distância entre os centros das par- 
tículas. 

O problema foi também estudado posteriormente 
por WOOD (5) e ROTENBERG 6), 

Um sistema deste tipo não pode ter realidade física, 
porque é um sistema a duas dimensões e porque o poten- 
cial intermolecular escolhido é extremamente simples. 
No entanto as conclusões a que se chegou são muito in- 
teressantes e merecem uma descrição pormenorizada. 

Na fig. 2 apresentam-se os resultados obtidos pelos 
referidos autores. As variáveis são a pressão reduzida 


pV 
J - (16) 
Nk T 
e o volume reduzido 
V 
- — — (17) 
V 


em que V, é a área duma configuração compacta hexa- 
gonal 


V —— No” (18) 
2 


A curva a cheio FV foi obtida por aplicação da 
teoria das células ou do volume livre («free-volumes), 
teoria que tenta traduzir as propriedades dos líquidos 
do ponto de visto da sua analogia com o estado sólido, 
o com a qual se obtém a equação 


t | [212 = - l (19) 
A figura 2 mostra uma boa concordância dos valo- 


res obtidos por Monte Carlo com a equação (19) para 
densidades elevadas (valores pequenos de 771). 
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Fig. 2 —- Resultados do método de Monte Carlo para a equação 
de estado dum sistema de discos rígidos: O o, wooD (5), com 
N=12;6 ! e WOOD !5), com N=48; +, METROPOLIS 14), 
Com N=224; À. ROTENBERGI8), com N= 224; as curvas 
a cheio representam a equação de estado segundo a teoria 
do volume livre (FV) e segundo a aproximação de REE e HOOVER 
à equação do virial P(3,3). 

(Reproduzido de !2) com a permissão da North Holland Pub. 
Co., Amesterdão). 


A curva a cheio P(3,3) representa a equação obtida 
a partir duma equação virial modificada, obtida por REE 
e HOOVER[7): 


es] + 1,8138077? (20) 


(1-0,3567803 + 7! + 0,021447 =?) 
(1-1,775171 : |! + 0,7787808: 


À equação do virial é aplicável ao estado gasoso, 
a a sua extrapolação para densidades da ordem dos 
líquidos corresponde a encarar o estado líquido do ponto 
de vista da sua analogia com o estado gasoso. Para den- 
sidades baixas (valores elevados de 771) pode observar- 
-se uma boa concordância entre a curva a cheio P(3,3) 
e os valores obtidos por Monte Carlo. 

O mais interessante na fig. 2 é, no entanto, a exis- 
tência de um intervalo no qual aparecem valores de pres- 
são reduzida mais elevada, que concordam razoávelmente 
com a equação virial, e outros valores de pressão redu- 
zida mais baixa, que se aproximam da equação proposta 
pela teoria das células. WOOD(2) designa este intervalo 
por intervalo de confusão e vem indicado por setas na 
fig. 2. Para melhor compreensão do que se passava 
nesta zona de densidade WOOD(2) estudou pormeno- 
rizadamente as experiências efectuadas. Cada determi- 
nação, constando de umas cadeias de n passos, foi divi- 
dida em M conjuntos de p passos cada, de tal modo que 
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M.p = n. Calcularam-se os valores médios para cada um 
dos conjuntos de passos da função de distribuição radial 
cumulativa, função relacionada com a função de distri- 
buição radial simples por uma equação integral. 

Na fig. 3 apresenta-se o resultado deste estudo, 
representando s o n.º de identificação de cada um dos 
M conjuntos referidos. 


0.50 o COESaS y E ES 
| 
| 
045: 
Es , 
pi " 
É 0s0/ | |A) RA | 
mo val |] + 
| 
| + 
035 , | | 
Rh”, 
Pr 
0.30 de — - — + 
(6) 20 40 60 Bo 
s 
Fig. 3 —- Carta de controle duma experiência com 48 discos 
rigidos, para - = 1,316 (5). Em ordenadas representa-se a 
função de distribuição radial cumulativa. Nesta experiência 


p= 19200 e M = 82. 
(Reproduzido de (2) com a permissão da North Holland Pub. 
Co. Amesterdão). 


Este gráfico, ou carta de controle, não tem clara- 
mente uma aparência aleatória mas, pelo contrário, apre- 
senta dois níveis de valores para a função representada, 
perfeitamente distintos. As médias obtidas sobre todo 
o conjunto de valores têm um significado dúbio e se- 
gue-se o processo de determinar os valores médios sepa- 
radamente para os valores nos dois níveis diferentes. 

Examinando um grande número de configurações, 
tais como as apresentadas na fig. 4, WOOD(2.5) con- 
seguiu estabelecer uma correlação entre o nível da carta 
de controle e a ocorrência de certos tipos de arranjos 
de partículas. 

Na fig. 4, a configuração s = 39, razoavelmente 
ordenada, corresponde a uma «pressão baixa» na fig. 2, 
enquanto as configurações correspondentes a uma «pres- 
são elevada», ou são desordenadas, como é o caso para 
s = 5, ou apresentam um arranjo curioso, de sete linhas 
de 6 discos e um buraco cada uma, 

Por outro lado, os deslocamentos das partículas em 
estados correspondentes aos «níveis baixos» são bastante 
pequenos em relação às distâncias aos vizinhos mais 
próximos, mantendo-se durante um n.º de passos elevado 
em volta das suas posições médias, correspondendo a 
uma estrutura hexagonal, enquanto nos estados corres- 
pondentes aos «níveis altos» se verificam grandes deslo- 
camentos. 

As transições entre o «nível baixo» e o «nível alto» 
e vice-versa são razoávelmente frequentes perto do ponto 
de densidade mais baixa do intervalo de confusão, mas 
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Fig. 4 —- Configurações com s = 5,7,39 da experiência apresen- 


tada na figura 3. O rectângulo interior é a célula unidade de 
volume V. As configurações com s= 5,7 estão associadas ao 
fluido no nivel mais alto, enquanto com s = 39 está associada ao 
nível mais baixo. 
(Reproduzido de (2) com a permissão da North Holland Pub. 
Co. Amesterdão). 


tornam-se rápidamente menos frequentes para densida- 
des crescentes, não se observando nenhuma para 7<1,3. 

O aparecimento deste intervalo de confusão num 
sistema de discos rígidos sugere a existência duma tran- 
sição de fase sólido-fluido nesta zona. 


4.2 — Esferas rigidas 


Se passarmos ao estudo de sistemas de esferas 
rígidas, cuja única diferença em relação ao de discos 
é o de se trabalhar em 3 dimensões, representando 
agora r, na expressão t15), o módulo dum vector tridi- 
mensional chegamos à conclusão que os resultados 
obtidos são semelhantes. 

Na fig. 5 apresentam-se os resultados obtidos por 
WOOD et ali8) ROSENBLUTH e ROSENBLUTH(9), 
ROTENBERG(S) e ainda os de ALDER e WAINWRI- 
GHT, (10) estes obtidos pelo método da dinâmica mo- 
lecular. 

Neste caso V, é o volume duma configuração com- 
pacta, cúbica de faces centradas, 


(21) 


À equação correspondente à curva a cheio P(3,3), 
obtida por REE e HOOVER,!7) é agora 


(22) 


1 + 0,063507 b: + 0,01732 (b:) 


4 — 0,561493 b: + 0,081313 (bo? 


Fig. 5 — Equação de estado de esferas rígida, segundo o mé- 
todo de Monte Carlo, As curvas a cheio têm significado seme- 
lhante ao da fig. 2; +, ROSENBLUTH e ROSENBLUTH (9), 
LI N=32. (4 N=-108. [E N- 256. ALDER e WAINWRIGHT 
(10), o N=32, 6 N=256, X N=32, WOOD et al (8), 
A N=32, 4 N- 864, ROTENBERG (8). 

(Reproduzido de (2) com a permissão da North Holland Pub, 
Co. Amesterdão). 


em que 


bs = 2,961922 71 (23) 


Novamente se verificam os fenómenos que se verifi- 
cavam para os sistemas de discos rígidos, cartas de 
controle com pontos corespondentes a «níveis altos» e 
pontos correspondentes a «niveis baixos», no intervalo 
de confusão, que aparece aqui na zona de volumes redu- 
zidos de 1,5 a 1,6, aproximadamente, 

HOOVER e REE (11) utilizando um modelo de célu- 
las, que são expandidas até densidades baixas, consegui- 
ram provar que a transição de fase sólido-liquido obser- 
vada para estes sistemas de esferas rígidas, é uma tran- 
sição de 1.º ordem, em que portanto a pressão se man- 
tém constante. 

Este resultado constitui talvez a contribuição mais 
importante dos estudos de simulação em computador 
para a compreensão das propriedades de equilíbrio do 
estado líquido. 

É claro que um modelo de potencial tão simplifi- 
cado como o de esferas rígidas, com uma parte infini- 
tamente repulsiva e sem parte atractiva, não pode con- 
duzir a uma descrição completa do comportamento dos 
líquidos, No entanto há provas experimentais de que 
a estrutura de um líquido, dada pela função de distribui- 
ção radial é práticamente independente da natureza das 
forças atractivas intermoleculares: a densidades corres- 
pondentes, a estrutura de sistemas metálicos e iónicos 
é quase a mesma do que a de sistemas de gases ra- 
ros (20), 
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Por outro lado, a ideia de que o comportamento 
ce um fluido denso é essencialmente determinado pela 
parte repulsiva do potencial intermolecular, sendo o po- 
tencial atractivo tratado como uma perturbação corres- 
pondente a uma energia coesiva média que mantém o 
sistema a uma densidade elevada, tem-se revelado muito 
importante na previsão teórica das propriedades de flui- 
dos densos. 


Várias teorias têm aparecido (ver ALDER e HOO- 
VER (20), BARKER e HENDERSON (12), RIGBY (21)) 
conduzindo a resultados em concordância muito razoável 
com os resultados experimentais para fluidos simples. 


O modelo físico que serve de base a estas teorias 
é provavelmente o mais conhecido modelo de fluido e 
foi imaginado há já mais dum século por Van der 
Waals (22), Ao escrever a equação de estado, p=f (V,7T), 
dum fluido, Van der Waals começou por separar os 
efeitos das forças repulsivas dos das forças atractivas. 


Uma vez que a lei dos gases perfeitos. pV = NkT, 
correspondente a um modelo físico em que não há inte- 
racções moleculares, Van der Waals introduziu dois tipos 
de correcções a esta lei: uma devida às forças repulsivas 
e outra às forças atractivas. Escrevendo a equação de 
Van der Waals na seguinte forma: 


NkT N 2 (24) 
ni = | =] 
V>-Nb V 
o primeiro termo à direita corresponde à correcção ce- 
vida às forças repulsivas, pela introdução de um parâme- 


tro b = - nº , correspondente à zona de influência 


duma esfera rigida de diâmetro 7. O segundo termo é a 
correcção devida às forças atractivas, partindo da hipó- 
tese de que, em fluidos densos, a resultante das forças 
atractivas exercidas numa molécula por todas as molé- 
culas vizinhas é nula, e portanto o potencial atractivo 
é constante. 


A principal objecção a pôr à equação de VDW é a 
de não prever a transição sólido-líquido. Embora Van der 
Waals procurasse uma equação que traduzisse o com- 
portamento dum fluido na transição líquido-gás, se o 
modelo físico fosse correcto seria razoável esperar que a 
equação previsse também a transição sólido-líquido. 


Os cálculos de Monte Carlo e da Dinâmica Mole- 
cular vieram demonstrar que não é o modelo físico que 
está errado, mas sim a aproximação feita à equação de 
estado de esferas rígidas que não é adequada. Na reali- 
dade se na equação (24) eliminarmos o termo — a(N/V)2 
deviamos obter a equação de estado de esferas rígidas. 
Ora a aproximação de VDW só concorda com os valores 
obtidos no computador para baixas densidades e, por- 
tanto, para densidades elevadas é necessário, para tornar 
a equação mais precisa, substituir o primeiro termo do 
membro direito da equação (24) por um termo adequado 
aos valores obtidos no computador para toda a gama 
de densidades. 


Houve várias tentativas para modificar este termo, 
como por exemplo, a de Frisch-Longuet-Higgins-Widom, 
a de Guggenheim e a de Scott 123), 
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Um dos resultados mais importantes da simulação 
em computador de sistemas de esferas rígidas foi o re- 
nascimento do modelo de Van der Waals, 


4.3 — Outros potenciais 


Das considerações anteriores conclui-se que o com- 
portamento dos líquidos, não pode ser descrito por um 
modelo de potencial tão simplificado como o de esferas 
rígidas. 

WOOD e PARKER 113) utilizaram o potencial de 
Lennard-Jones (6-12). 


G 12 G 6 (25) 
ul) —- 4: Ra e ) 


A forma aproximada deste potencial é representada 
na fig. 6 


U(r) 


Fig. 6 — Potencial de Lennard-Jones (6-12). 


O artigo de WOOD e PARKER constitui a primeira 
aplicação de método de Monte Carlo ao estudo dum 
fluido real. 

VERLET e LEVESQUE (14) e Mc DONALD e SIN- 
GER(15) estudaram também as propriedades dum siste- 
ma com um potencial intermolecular deste tipo e os 
seus resultados concordam bem com os valores das 
propriedades do àArgon, pelo menos para temperaturas 
entre o ponto triplo e ponto crítico. 

LEVESQUE e VIEILLARD-BARON (16) efectuaram 
cálculos para vários tipos de potencial e chegaram à 
conclusão que o potencial de Lennard-Jones (6-12) é o 
que se adapta melhor às propriedades de Argon. 

A escolha do Ar para termo de comparação justi- 
fica-se pelo facto das suas moléculas serem monoató- 
micas e esféricas. Deste modo o potencial de interacção 
molecular não contém certamente nenhum termo depen- 
dente da orientação recíproca das moléculas e deve 
ser função unicamente da distância entre elas, Por outro 
lado, como é uma molécula relativamente pesada, são 
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desprezáveis os efeitos quânticos que conduzem a des- 
vios em relação às leis da Termodinâmica Estatística 
clássica, 

A concordância com os valores experimentais 
das propriedades do Ar constitui portanto um teste muito 
efectivo da validade dum determinado modelo de po- 
tencial. 

No entanto, uma análise cuidadosa dos dados ex- 
perimentais mostra que o verdadeiro potencial entre duas 
moléculas de Ar tem um fosso de potencial mais pro- 
fundo e mais abrupto que o do Lennard-Jones (6-12). 
Para baixas densidades (gás pouco denso) o potencial 
que está mais de acordo com as propriedades do Ar 
é o potencial de Kihara, em que se admite a existência 
dum núcleo rígido de diâmetro a: 


q—-—a 12 — a. 6 


) — ) r>Zs 
f—:B Fr | 


IH 
s 


u (r) 


(26) 


II 
R 


u (r) r< 28 
Este paradoxo pode ser resolvido se admitirmos que a 
contribuição de potenciais entre N corpos, não aditivos, 
tem para densidades da ordem das do Ar líquido, uma 
contribuição não desprezável. 

Me DONALD e SINGER |3) efectuaram cálculos com 
o potencial de Kihara, fizeram estimativas, também pelo 
método de Monte Carlo, das forças entre 3 dipolos 
e, somando os resultados, obtiveram valores em muito 
melhor concordância com os valores experimentais para 
o Ar que os obtidos só pelo potencial de Kihara, embora 
a concordância não seja ainda boa como a obtida com 
o potencial de Lennard-Jones (6-12). 

Este funciona portanto como um potencial efectivo 
entre duas moléculas, 


4.4 — Simulação de moléculas lineares em cadeia 


Os métodos de Monte Carlo têm sido também apli- 
cados nos últimos anos ao estudo de moléculas polimé- 
ricas lineares. 

Os primeiros estudos foram realizados por WALL 
et al (24), para cadeias de segmentos com um potencial 
de interacção de esferas rigidas. 

No caso de os segmentos interactuarem a distância, 
a energia configuracional é diferente de zero pelo que 
se aplica a expressão (1), uma amostra aleatória- 
mente escolhida entre configurações possíveis. No en- 
tanto, as dimensões da amostra têm que ser relativa- 
mente grandes para a convergência da média ser satis- 
fatória. Este foi o processo empregado por WALL e 
MAZUR (25) e outros 128,27) que obtiveram um relativo 
sucesso para potenciais de fosso rectangular. 

O uso do potencial de Lennard-Jones (6-12) está 
limitado à simulação de cadeias muito pequenas!28), 
pois envolve cálculos muito demorados. 

LAL(29), aplicou o método de METROPOLIS et al!4) 
a uma molécula linear em cadeia (polimero). 

O método utilizado foi o descrito na secção 2, mas 
convém ilustrá-lo para uma cadeia bidimensional numa 
rede hexagonal, como se mostra na fig. 7. 
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Fig. 7 — Simulação duma cadeia bidimensional numa rede 


hexagonal, 


Sistematizando temos: 

1) Escolhe-se um segmento ao acaso m, e consi- 
dera-se a nova configuração como reflexo na a linha defi- 
nida por m, da parte CB. 

2) Se a nova configuração ACB' conduzir a uma 
ocupação múltipla das células da rede escolhida, rejei- 
ta-se esta configuração e a nova configuração é ACB. 
Caso não conduza, aceita-se a nova configuração. 

Como alternativa deste processo temos o cálculo 
da diferença de energia configuracional de ACB e ACB., 
AU. Se AU - 0, a cadeia toma a nova configuração e se 
AU < O compara-se exp (A U/kT) com um número 
aleatório:, 0< E<'1, e se exp (AU/kT) —: aceita-se a 
nova configuração; caso contrário rejeita-se ACB', re- 
toma-se ACB, escolhe-se outro segmento p e assim su- 
cessivamente. Podemos esquematizar este último pro- 
cedimento num diagrama de blocos, fig. 8. 


Fig. 8 — Diagrama de blocos exemplificativo da simulação 
duma cadeia. 


Nos cálculos efectuados LAL admitiu a aditividade 
da interacção intersegmental e o potencial de LJ (6-12), 
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com :k= 88 K e r' = 3.34, determinados a partir do 
princípio de estados correspondentes, segundo o pro- 
cesso de PRIGOGINE /30), 

Foram suficientes 100 000 configurações para a 
convergência satisfatória das propriedades configuracio- 
nais médias duma cadeia de 30 segmentos, gerados 
aproximadamente em 1/2 hora num computador IBM 
360/50. 

A incerteza é da ordem de + 1,0% no cálculo de 
<UN>IN é o número de segmentos), e a forma das 
curvas <U, -/N em função de kT/: prevê transições de 
ordem superior nas cadeias. 

LAL(31) estudou também o volume excluído das 
cadeias, isto é, o volume não ocupado pelas cadeias 
devido à repulsão intersegmental, À convergência obtida 
nos cálculos foi óptima, acima de 200 000 configurações. 
Para uma sequência de 500 000 configurações um IBM 
360/50 levou aproximadamente 4 horas. 

Os estudos de Monte Carlo, devido a serem semi- 
-experimentais, podem utilizar potenciais mais realistas, 
baseados em dados experimentais da geometria e do 
comportamento energético das moléculas do polímero. 

Mais recentemente LAL (32), considerou um modelo 
mais realista para as moléculas de n-alcanos; entrou 
com um potencial rotacional das ligações C-C, com 
interações de curto alcance e de longo alcance. 

As interacções de curto alcance foram determinadas 
a partir dum potencial de Buckingham (exp-6), com os 
parâmetros de ABE et al (33). As de longo alcance fo- 
ram calculadas admitindo um potencial de LJ, com 


:/k=88,4K e r' = 4,19 A, valores estes derivados do 
calor latente de sublimação e das densidades a OK do 
polietileno. 

Foram estudados o comportamento conformacional 
das ligações a várias temperaturas e a capacidade calo- 
rífica, esta para pôr em relevo as transições já observa- 
das só com o potencial único 129), 


5 — CONCLUSÃO 


Nos últimos anos, o desenvolvimento da capacidade 
de cálculo e memória dos computadores tem permitido 
a realização de cálculos para sistemas cada vez mais 
complicados e com potenciais de interacção cada vez 
mais complexos. 

São de salientar o trabalho de Mc DONALD [17,34] 
sobre misturas de líquidos simples, o de RAHMAN e 
STILLINGER (18) para o estudo das propriedades da 
água utilizando o método da dinâmica molecular e os 
cálculos de BARKER, FISHER e WATTS (19) sobre as 
propriedades do argon, 

Quanto às cadeias poliméricas, o desenvolvimento 
também tem sido grande, e será futuramente aplicado 
a interacções molécula-superfície, estudo de transições 
de ordem superior e de 1.º ordem. 

O método de Monte Carlo, devido ao seu carácter 
semi-experimenta!l promete ainda maior desenvolvimento, 
pois os resultados obtidos são, na maior parte das vezes, 
melhores do que os resultados obtidos por via teórica. 
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O Plano Geral de Actividades de Formação e Aperfeicoamento do INII para 1972 - 73, acaba 
de ser divulgado pela Divisão de Formação e Aperfeicoamento. 

A programação das acções assenta em quatro programas de base: Estratégia e Estrutura 
das Organizações — vai realizar-se de 6 a 10 de Março; Estratégia Geral da Empresa e Exercício 
de Direcção — terá início a 19 de Fevereiro; O Curso Superior de Management, iniciando-se a 15 
de Janeiro e a Direcção e Gestão Modernas da Pequena e Média Empresa — será iniciada a 26 
de Janeiro. 

No Plano Geral seguem-se os Ciclos de Desenvolvimento em Direcção; Ciclo de Marketing, 
constituído por 4 módulos obrigatórios e 5 especializados. O primeiro terá início a 19 de Março 
— Ciclo Direcção Financeira, formado por 3 módulos obrigatórios e 2 especializados — o pri- 
meiro inicia-se a 19 de Março — Ciclo Direcção de Pessoal — 2 módulos obrigatórios e 2 espe- 
cializados — início a 12 de Março — Ciclo Management e Informática — formado por 3 seminá 
rios a comecar dia 23 de Abril e Ciclo de Análises e Decisão — 3 seminários, com início a 2 
de Abril. 

A completar as actividades para 1972 - 73, o Plano Geral refere ainda alguns Programas 
Específicos em colaboração com a Direcção de Servicos da Reforma Administrativa, com a 
Direcção-Geral dos Combustíveis e Banco de Fomento Nacional, 6 Sessões de Introdução e 9 
Seminários Especializados. 

Vão ainda realizar-se um Simpósio sobre «A Formação e o Desenvolvimento dos Recursos 
Humanos» em Abril e as Acções no Ultramar em Luanda e Lourenço Marques, em colaboração 
com as Associações Industriais de Angola e Moçambique. 

Qualquer informação poderá ser solicitada a —INII — Divisão de Formação e Aperfeicoa- 
mento — Rua Filipe Folque, n.º 2-6.º Dto. — Lisboa -1 — Telefs. 53 36 24/5/6. 
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J. J. DELGADO DOMINGOS Cc. D. U. 600 
A CRISE DO AMBIENTE 
Técnica N.º 417 — XLVIIL — 11-1972, p. 57-64 


Abordam-se discutem-se problemas ligados à crise do am- 
biente e clarificam-se aspectos e afirmações correntes a 
partir de princípios fundamentais da Física. As possibili- 
dades e limitações da Tecnologia são perspectivadas. Algu- 
mas implicações sociais, económicas e políticas são in- 
terpretadas mostrando-se como as questões de evolução e 
sobrevivência dos povos são marcadamente diferentes con- 
soante o seu poder e grau de desenvolvimento económico. 


C. D. U. 551.311.22/620.199/691.2 : 5392.15 
539.215 
L. AIRES - BARROS 


ANALISE LABORATORIAL DA INFLUÊNCIA DA GRANULA- 
RIDADE NA ALTERABILIDADE DAS ROCHAS 
Técnica N.º 417 — XLVIIL — 11-1972, p. 65-70 
Faz-se o estudo da influência da granularidade das rochas 
na sua alterabilidade. Para isso utilizam-se extractores de 
Soxhlet, tendo-se ensaiado três tipos texturais de uma 
mesma rocha ígnea (sienito nefelínico). A partir da mobi- 
lidade quina dos principais catiões constituintes dos mi- 
nerais das rochas e das perdas de peso de cada tipo de 
material ensaiado, procura-se adaptar um Índice de alte- 
rabilidade. Tenta-se avaliar o peso a atribuir à textura de 
uma rocha ao se usar uma expressão geral da alterabili- 

dade das rochas. 


C. D. U. 620.17 : 624.073,12 
MANUEL AMÉRICO GONÇALVES DA SILVA 
RESPOSTA ELASTICA DUMA LAJE CIRCULAR 
Técnica N.º 417 — XLVI| — 11-1972, p. 71-77 


Investiga-se a resposta elástica duma laje circular de bor- 
dos encastrados, sujeita a uma carga transversal distri- 
buída, utilizando diversas teorias que se comparam. Ana- 
lisam-se vários modos de aproximar o apoio encastrado 
e estuda-se a influência da espessura relativa da laje 
na acuidade dos resultados. 


J. F. BORGES DA SILVA C. D. U, 621.3.011.1 
ALGORITMO PARA A ANALISE TOPOLÓGICA DE REDES 
Técnica N.º 417 — XLVII — 11-1972, p. 79 - 86 


São expostos os fundamentos teóricos e a organização ló- 
gica dum algoritmo destinado à análise automática de redes 
pelo método topológico, usando um computador digital. O 
algoritmo pode ser usado indiferentemente com informação 
de incidência relativa a circuitos ou a cortes, sendo pos- 
sivel a análise de redes com ligação magnética e incluindo 
elementos não recíprocos activos ou passivos. 


ACÁCIO DA CONCEIÇÃO C. D. U. 691.327 : 620.197.14 


PROTECÇÃO DE DEPÓSITOS E OUTRAS ESTRUTURAS 
EM BETÃO POR MEIO DE REVESTIMENTO 


Técnica N.º 417 — XLVIL — 11-1972, p. 87 - 96 


Devido ao uso cada vez maior do betão de cimento 
portland tanto no domínio da construção e obras públicas 
como no das técnicas e operações industriais, incluindo a 
armazenagem de produtos químicos e alimentares os mais 
diversos, estão continuamente surgindo novos problemas 
do seu ataque químico, ao mesmo tempo que vão apare- 
cendo no mercado novos produtos para a sua protecção 
cuja adaptação a dadas condições de corrosão deve ser 
antecipadamente confirmada em cada caso prático espe- 
cífico. 

No presente trabalho bibliográfico fazem-se algumas con- 
siderações gerais sobre o ataque químico do cimento por 
produtos de uso industrial e afloram-se os processos mais 
viáveis para proteger as superfícies de betão por meio de 
revestimentos anticorrosivos. No capítulo final resume-se 
a matéria delineada e indica-se a principal bibliografia a 
que, em princípio, se poderá recorrer para definir e esta- 
belecer o processo, ou processos mais aconselháveis, para 
revestir o betão de cimento portland em face de cada 
produto industrial ou alimentar. 


JOAO REIS SIMÕES C. D. U. 512.8/9 
FORMAS MATRICIAIS NO CALCULO DE AFINORES (ll) 
Técnica N.º 417 — XLVIIL — 11-1972, p. 97-99 


Adaptam-se a algumas operações sobre afinores, as for- 
mas matriciais. 


C. D. U., 54.02 - 14: 681.117.8: 519.2 


C. A. NIETO DE CASTRO 
M. NUNES DA PONTE 
V. A. MEIRA SOARES 


APLICAÇÕES DO METODO DE MONTE-CARLO 
AO ESTUDO DE FLUIDOS DENSOS 


Técnica N.º 417 — XLVIL — 11-1972, p. 101 - 110 


Apresentam-se em revisão algumas aplicações do método 
de Monte Carlo ao estudo de fluídos densos, abordan- 
do-se alguns dos modelos utilizados na simulação das 
propriedades dos fluídos densos, como os discos rígidos, 
esferas rígidas e Lennard - Jones; resume-se ainda a apli- 
cação à simulação de cadeias poliméricas lineares. 
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ACÁCIO DA CONCEIÇÃO  U. D.C. 691.327: 620.197.14 


PROTECTION OF STORES AND OTHER REINFORCED 
CONCRETE STRUCTURES BY MEANS OF COATINGS 


Técnica N.º 417 — XLVII — 11-1972, p. 87 - 96 


Owing to the increasing use of portland cement concrete 
not only in buildings and public works but in process engi- 
neering operations -—- including the storage of chemical 
and food products —, new problems are constantly arising 
in connexion with chemical attacks on these structures. 
At the same time new products for protection against these 
attack are constantly appearing in the market and their effi- 
ciency in given corrosion problems must be confirmed 
beforehand in each particular case. 

The present bibliographical report begins with some gene- 
ral considerations on the chemical attack of cement by 
industrial products, and briefly introduces the most feasible 
methods of protecting concrete surfaces by means of cor- 
rosion-resistant coatings. The final chapter sums up the 
subject of he report, and lists main references which in 
principle should be used for choosing the most suitable 
coating procedures for protection of portland-cement con- 
crete against given industrial and food products. 
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with 
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are adapted to some operations 
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APPLICATIONS OF MONTE-CARLO METHOD TO THE 
STUDY OF DENSE FLUIDS 


Técnica N.º 417 — XLVIIL — 11-1972, p. 101 - 110 


A review of some applications of the Monte Carlo methods 
to the study of dense fluids is presented. Some of the 
models used in the simulation of the dense fluid properties 
are discussed, such as the hard rod, the rigid sphere and 
the Lennard-Jones models. The application to the simulation 
of the chain molecules is presented as well. 


TÉCNICA 417 


J. J. DELGADO DOMINGOS 
THE ENVIRONMENTAL CRISIS 
Técnica N.º 417 — XLVIIL — 11-1972, p. 57 - 64 


U. D. €. 600 


Some problems of the environmental crisis are presented 
and discussed and current statements clarified from fun- 
damental physical principles. Possibilities and limitations 
of modern technology are discussed. Political, social and 
economic implications are analysed and it is shown how 
strongly are the survival possibilities of people connected 
with its power and economic development. 
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LABORATORY TEST ON THE INFLUENCE OF THE GRA- 
NULARITY ON ROCK ALTERABILITY 


Técnica N.º 417 — XLVII — 11-1972, p. 65- 70 


The study on the influence of the granularity of rocks on 
their alterability is intented. 

Laboratory tests were conducted by means of Soxhlet 
extractors using three different grain-size types of the 
same rock (nepheline syenite). 

Using the chemical mobility of main cation of the mineral 
constituents and also the weight lost for each rock type 
analysed, and adaptation to an alterability index is inten- 
ted. This study attemps to know the textural factor of 
a rock which may be employed in a general equation of 
rock alterability. 
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ELASTIC RESPONSE OF A CIRCULAR SLAB 
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One investigates the elastic response of a circular slab 
with fixed edges, subject to a transverselly distributed load, 
using different theories which should he tested. One ana- 
lyses different ways of approximating the fixed support and 
one studies the influence of the relative thickness of the 
siab in the accuracy of the results. 
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The theory and logical structure of an algorithm for auto- 
matic network analysis by digital computer using topolo- 
gical methods are described, The algorithm may be used 
indifferently with incidence data given for circuits or cut- 
-sets. The analysis of networks exhibiting magnetic coupling 
and including active or passive non-reciprocal elements is 
possible. 
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— F. P. Kotchnev 

Rail Int. 3 (2) 1972, 102 - 113 


C. D. U. 663.6 


Nuevo Método en la depuración de aguas residuales — 
Justo Leácer 
informes de la Construcción, 24 (238), págs. 69 - 72, 
Março 1972 


C. D. U, 666.912: 621.316.7-523.8 : 681,3 
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Seabra 
Memória 388 do LNEC, 1971 
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C. D. U. 669.2.886 : 669.5 
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As instalações eléctricas em Portugal — José Ribeiro 
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Empresa de Sondagens e Fundações 
TEIXEIRA DUARTE, L.” 


o même LISBOA 


SONDAGENS GEOLÓGICAS 
CAPTAÇÃO DE ÁGUAS SUBTERRÂNEAS 


CONSOLIDAÇÃO DO SOLO 
INJECÇÕES DE CIMENTO 


FUNDAÇÕES DE TODOS OS TIPOS 
CAVES E TUNEIS 
BARRAGENS E PORTOS 


(UM QUARTO DE SÉCULO DE ESPECIALIZAÇÃO TÉCNICA) 


6170 


Disjuntor 


Protege pessoas e animais 


m contra tensões de contacto excessiva- 
mente elevadas, em aparelhagem 
deficientemente isolada 


m em casos de troca do condutor 
principal pelo condutor de protecção, 
ou ainda, troca entre o condutor de 
protecção e o neutro 


m duma maneira geral, portanto, sempre 
que uma corrente perigosa passa, 
através do corpo de pessoas ou de 
animais, à terra 


Siemens 


Garante 


E grande protecção humana 


ms em caso de terra fraca, uma supersen- 
sibilidade (dispara para correntes 
de defeito entre 15 e 30 mA segundo 
prescrições VDE 0100 e dentro de 
30 ms) 


m elevado poder de corte 


m longa duração, cerca de 20.000 mano- 
bras à corrente nominal 


O 


“a, 
Casa 


A, 


ça 


SIEMENS 


melhor sistema 
de proteccao humana 


de corrente de defeito 


Aplica-se 


como interruptor geral em habitações, 
hospitais, oficinas, fábricas, instalações 
provisórias, obras e em todos os casos em 
que haja elevado grau de humidade ou 
deficientes condições de isolamento. 


Exponha os seus problemas especiais à 
secção técnica da 

Siemens-Companhia de Electricidade, SARL 
Av. Almirante Reis, 65-tel. 526921-Lisboa 1 
R. das Carmelitas, 26-2.º - Tel. 28943 - Porto 
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